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Abstract 

Generalizing the mapping between the Potts model with nearest neighbour interaction and the six 
vertex model, we build a family of "fused Potts models" related to the spin fc/2 U 9 su(2) invariant 
vertex model and quantum spin chain. These Potts models have still variables taking values 1,...,Q 
{VQ = 9 + Q 1 ) but they have a set of complicated multi spin interactions. The general technique to 
compute these interactions, the resulting lattice geometry, symmetries, and the detailed examples of 
k — 2,3 are given. 

For Q > 4 spontaneous magnetizations are computed on the integrable first order phase transition 
line, generalizing Baxter's results for k — 1. 

For Q < 4, we discuss the full phase diagram of the spin one (k = 2) anisotropic and U g su(2) invariant 
quantum spin chain (it reduces in the limit Q = 4 (g = 1) to the much studied phase diagram of the 
isotropic spin one quantum spin chain). Several critical lines and massless phases are exhibited. The 
appropriate generalization of the Valence Bond State method of Affleck et al. is worked out. 
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1 Introduction 



Numerous families of integrable lattice models have recently been exhibited. In general, these families 
obey the following pattern. One starts with an algebra like sl(n) and one of its representation say p. Using 
quantum group technology a solution of the Yang Baxter equation acting in p ® p can be found, which 
encodes the Boltzmann weights of a vertex model. These weights are trigonometric and they depend on 
the quantum group deformation parameter q and the spectral parameter. The degrees of freedom are the 
weights of the representation p, and they sit on edges of usually the square lattice. Using the quantum group 
symmetry the model can also be reformulated as a solid on solid model whose degrees of freedom are 
highest weights and sit now on the vertices of usually the square lattice. For q a root of unity, the model 
truncates and a restricted sos model can be defined. 

In the particular case of su(2) spin 1/2, there exists, besides the six vertex and sos model, a third 
"equivalent" model, the Potts model Qj. Its existence and relations with the first two have a precise 
meaning in terms of Temperley Lieb algebra representation theory. The purpose of this paper is to show 
that for higher spin, there is also a Potts model naturally associated with the vertex and sos models. This 
model still uses the same set of variables a = 1, . . . , Q but involves interactions that are more complicated 
than the nearest neighbour coupling of the spin 1/2 case. That such models exist has been known in principle 
for a long time fll9[ but their precise definition, and the general algebraic formalism to build them, are new 
to the best of our knowledge. While we were working on that construction, we became aware of the work 
of Nienhuis Q where the first member of the hierarchy, a Potts model with nearest and next to nearest 
neighbors interactions related to spin 1, has already been exhibited. However the techniques used in || and 
in our work are very different. 

The construction of what we shall call Potts model is discussed in details in section 3, after ele- 
mentary reminders about k = 1 are given in section 2. From another perspective, it gives an interesting 
light to the differences between integer and half integer vertex models or quantum spin chains (we can in 
particular reformulate the spin k/2 su(2) spin chains problems in terms of Q — 4 states Potts models). For 
k greater than one, only submanifolds of the full parameter space are integrable. In section 4 we discuss the 
simplest integrable line for Q > 4 where a first order phase transition takes place. We compute spontaneous 
magnetizations by generalizing Baxter's calculation for spin 1/2. 

We discuss in section 5 the full phase diagram of the k = 2 model, mainly in terms of the corresponding 
one dimensional quantum spin chain. This phase diagram restricts in the case Q = 4 to the widely studied one 
for the spin 1 su(2) hamiltonian. Critical properties are studied with particular emphasis on the integrable 
lines. We also construct the valence bond states for the P2 projector for arbitrary Q thereby extending 
the construction in JtJ to anisotropic systems (with quantum group symmetry). In opposition to [Q, phase 
transitions are encountered as Q varies. 

Two appendices are included. Appendix A contains the explicit expression for the Boltzmann weights of 
the k = 3 Potts model. Various loop model interpretations for the family of Potts model, in particular for 
k = 2, are reviewed briefly in appendix B. 
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2 The Q-state Potts model and staggered six-vertex model 



We begin by reviewing the Q-state Potts model on a square lattice. We consider a cylinder made of an 
/ x t rectangular strip C with free boundary condition on the top and bottom rows, and periodic in the time 
direction. The partition function of the model with horizontal and vertical couplings K\ and K 2 is given 

by§ 

z Poaa = j2U eKlS '" i+KaS ' i ' i > t 2 - 1 ) 

Wi} (ij) 

where the product is over all neighboring horizontal and vertical links < ij > and <7i assumes values 1 to Q. 

The column to column transfer matrix rp tts can be expressed as a product of local transfer matrices 
X 2 {-i and X 2 i which add, respectively, a horizontal and a vertical link to the lattice. We have 



i=l i=l 

and 



rpotts = Q W2 J] X 2i-! 1] X 2i (2.2) 
i=l 

^Potts = tr(r) 4 (2.3) 



with 

X 2 i = 1 + x 2 e 2 t , 

where we define 

e K k -l 



(2.4) 




Xk = -qTJT for fe=1 ' 2 

and the operators e^'s which propagate {ci} to {cr^} in the time direction, have matrix elements given 

(2.5) 

= Q 1 ^ 2 Y\j^Li5 ajCy ' j 1 
They satisfy the following relations (dropping the subscripts from now on) 

= ei, (2.6) 
= ;|i-i|>2. 

The algebra generated by them is known as the Temperley Lieb algebra. In the Potts representation the 
following trace properties hold|j| 

tr[w(l,ei,...,e i _i)e l ] = Q _1 / 2 tr w(l, e\, . . . , e,_i) , 

tr(ei) = Q-^Hrl, and (2.7) 
trl = Q l , 
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where w is any word in 1, ei, . . . , The Potts model exhibits a duality transformation implemented by 

rewriting the local transfer matrices as follows 



X2i-1 
X 2 % 

and interchanging the roles of e2i-\ and e 2 i, 

( e 2i) CT)CT < 

(e2i-i) CT v 



sci(l + x x 1 e 2 i-i) 
+ e 2l ) , 



n 1 / 2 s U l S ' 

V u a-i<Ti + llij=l u cr j a-' j , 



(2.8) 



(2.9) 



which amounts to redefining X 2 %— i and X2i to be the local operators that add a vertical and a horizontal 
link to the Potts lattice respectively. This alternative interpretation of the roles of and e 2 i does not 

alter the algebraic relations satisfied by them . We thus have 

Zpott s (:r 1,2:2) = [xxx 2 ) lt Z-p^^x^ 1 ,X^ X ) (2.10) 

where the prime denotes the dual lattice C which differs from C by boundary effects only, (see fig.(l)) 
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Figure(l) Geometry of the lattice and its dual 



The duality transformation 



xi< — > Xz (2.11) 

relates high temperature to low temperature phase. Ignoring the difference in the boundary, the model is 
self-dual at 

x lX2 = 1 , (2.12) 
and by standard argument |Q this is a line of phase transition (first order fo Q > 4, and second order for 
Q < 4). 

It is well known that the model can be mapped to the six-vertex ( referred to as T\ here ) model by 
assigning arrows on the surrounding polygons of the clusters formed by Potts variables <7j that have the 
same colors |l], |j| in the high temperature expansion. The procedure is more transparent from an algebraic 
point of view. Consider the tensor product of 21 copies of spin-i representations of U 9 sl(2), then can be 
represented in this vector space as 



( 



V 







\ 

- 1 -1 

-1 q 

/ 



e C 6 



(2.13) 
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where the matrix is proportional to the spin-0 projector of the i th and (i + l) th spin-i representation, it is 
easy to verify that the above indeed satisfies ( O ) with 



and the trace properties ( 2.7) can be reproduced by introducing the Markov trace defined as 



where 



tr 



w(l,ei 




tr ( q 2sZ a, 



tr I q 



(q + q- 1 ) 

(q + q- 1 ) 



tr 



tr [q 



q 2S w(l,ei, 

2S Z ' 



, e»_i) 



(2.14) 



and 



-3 



1 . 



The operator Xi, with e, given as ( 2.13| ), defines the vertex interaction of the staggered six- vertex model, 
that is its matrix elements encode the Boltzmann weights of the vertices with various colors for incoming 
and outgoing lines. The partition function is then given by 



= tr 



q 'vertex 



(2.15) 
(2.16) 



i=l 



with e<'s in the above defined by ( 2.13 ). Because of the same trace properties of the two representations, 
■^vertex and Zp tts are equal for integer Q. However, the former is defined as well for Q real (and coincides 
with the geometrical definition based on high temperature expansion). The Temperley Lieb generators 
commute with U 9 sl(2), so the six- vertex model has U 9 sl(2) symmetry. The mapping between vertex and 
Potts models is not one to one when q is a root of unity due to the boundary operator q 2S in the trace and 
the quantum group symmetry. It is actually one to one between the Potts model and only the subset of type 
II representations of U g sl(2) provided by the six- vertex model pH, |Iiq ]. 



3 Potts model formulation of the Tt vertex model 

3.1 The fusion procedure and mapping to Potts model 

We shall construct a family of Potts models which are related to the U g su(2) invariant vertex models 
based on spin-|. We call the latter vertex models. This name is non standard, and we do not know 
any more standard name to use instead (in M] Tk refers to the number of allowed vertices, 6 for k = 1, 
19 for k = 2, . . .). For the moment we therefore decide only about the symmetry, not the particular set 
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of interactions. The construction uses then ideas of the "fusion procedure" || (properly generalized) to 
reexpress each U q su(2) spin-fc/2 in terms of k copies of spin-1/2. A pairs of such 1/2 spins interact at 
vertices of the six- vertex model (T\ vertex model). Each such vertex is then translated into a Potts model 
interaction using the results of section 2. This provides finally a Potts model with complicated inhomogeneous 
interactions which we call the Potts model. 

This construction is best explained by explicit computation. First, the Boltzmann weights of a particular 
vertex with two incoming and two outgoing legs carrying U g su(2) spin-fc/2 variables are encoded in a matrix 
which we call for simplicity the vertex as well. It is an operator that acts on C fc+1 (g> C k+1 . The 
corresponding spin-| irreducible representations can be obtained from the g-symmetric tensor product of k 
copies of the spin-^ one. This is conveniently done using the g-symmetrizer defined by[[l4| 

q k(k-l)/2 



[k] q 

where 



^Ssn-E*" 1 '-'!!-^ (3 - L1) 



s(i) = q x l - e, . 

In the above formula, I a denotes the collections of indices in the nonreducible decomposition of the symmetric 
group element a in terms of transposition of neighbors t^+i, \I a \ indicates the cardinality of the collections 
and we also introduced the g-factorial where 



\ g \ = (n),(n -!),...(!), (3.1.2) 



g" - q~" 

and (n)n = t- 

q-q 



The symmetrizer can alternatively be written recursively as^j, ^] 

(*-!)« 



S k = S k -i 1 ttt — e k -i S k -i , (3.1.3) 

V ( k )q J 

which expresses S k solely in terms of products of Ti vertices. 
The Tfc vertex written in terms of Ti reads then 

SkSk[rk{ui)r k +i{u 2 ) ■ ■ ■ r 2k -i(uk)][rk-i(uk+i) ■ ■ ■ r 2 k-2(u 2 k)] ■ • • [ri{u k 2 _ k+1 ) ■ ■ ■ r k (u k 2)]S k S k , (3.1.4) 

where the first and last two S k s act respectively on the two in- and out- states, and we encoded the T\ 
weights in the matrix 

sin ti, 

r i (u j ) = l + — (3.1.5) 

sm(7 — Uj) 

with q = e %1 . Note that Tj(uj) is identical with the local operator Xi introduced in the previous section 



The construction ( 3.1.4 ) is illustrated graphically in fig. (2), which shows that the interaction between 



the two spin--| states is replaced by the interactions r k {u\) 1 ■ ■ ■ 7 r k (u k 2) between 2k spin-i states. It is not 



difficult to see that the operator ( 3.1.4 ) has U g su(2) symmetry since each of the Ti vertices regarded as an 
operator in C k+1 (g> C k+1 also has U g su(2) symmetry. This together with the fact that S k 's project onto the 
spin- | irreducible representation, implies that the operator ( |3.1.4 ) is indeed a T k vertex . It can therefore 
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be written as a linear combination of the projectors Pj ,j = 0, 1, • • • , k. Conversely it can be shown ( see 
also appendix B) that any Tk vertex weight can be written as the above using a particular set of m's. The 
number of independent parameters for the Tk vertex is equal to k (we factored out an irrelevant overall 
scale), which is much less than the number of itj's. But these are the most convenient parameters. 




The fused vertex 



The symmetrizers S k 



Figure(2) The fused i\ vertex; The parameters ui, ■ • ■ ,u k 2, etc appearing in the figure 

are associated with the Ti vertices that make up the Tj. vertex. The four S^'s that 
act on the in- and out- states arc shown on the rhs. 



Next, having written the Tk vertex in terms of Ti vertices and thus the Temperley-Lieb generators, we 
can use the various realizations of the Temperley-Lieb algebra. In particular, we shall consider the Potts 
model realization, which produces the Tk Potts model we set out to construct. The Potts model lattice 
correspondi ng to the Tk vertex is built up as follow s: To each Ti vertex in ( 3.1.4 ), we substitute the 
expression ( |2.5|) or ( 2.E ), and call the operator ( 3.1.4 ) in this representation W(u), and the matrix 



elements of W{u) induced by the in- and out- states of the vertex will become the Boltzmann weights of 
the Tk Potts model. Graphically we assign either a horizontal or vertical link to each spin-i vertex rj as 
shown in fig. (3), W(u) therefore corresponds to an operator that acts on k or k + 1 Potts variables, and is 
represented graphically as a connected collection of horizontal and vertical links, which shall be referred to as 
the fundamental block Qj~ hereafter. The lattice C is then constructed by replacing all the vertices by 
these fundamental blocks, we shall elaborate on this point in next few sections. Note that there are exactly 
two possible mappings to the Potts model links which originate from the two possible choices of assigning 
links to the generator, namely ( 2.5 ) and ( 2.£ ) . When ( 2.5 ) is used, a horizontal ( vertical ) link is associated 
to rj with odd (even) subscript, whereas when (2.9) is used, the roles of vertical and horizontal links are 
reversed. We shall denote these two choices of mapping as convention A and B respectively. It will be 
shown that they are related by duality transformation which is inherited from that of the Ti Potts model. 
The relation between the Tk vertex and Potts models generalizes that of the k = 1 case. In particular, 



G 



for the vertex model whose lattice has the geometry of a cylinder with periodic time boundary condition, 
the partition function defined as in ( |2.15| ) is equal to that of the corresponding Potts model when Q is an 
integer. The Markov trace in ( [2.15| ) again restricts the domain of the mapping to type II representations of 
U g su(2) provided by the vertex model. 




C '' d 

Figure(3) Mapping of the Y\ vertices to Potts model links. In the first figure on the rhs, 

Tn— 1 ixii) i s replaced by the horizontal ( vertical ) Potts model links that connects 
Oi and a i ( a i and o" 4+1 ), the dual Potts model is given by the second figure on 
the rhs where the above vertex is mapped to a vertical ( horizontal ) links 
connecting Ui and Gi+x ( cr i+ i and o" i+1 ) 

3.2 The local interaction and its dual 



We shall first work out explicitly the weight W(u) for k = 2 and 3, and construct graphically the Potts 
model representation of W(u), ie. the fundamental blocks Qi and Q 3 . The construction of the fundamental 
block is then generalized to arbitrary k. 

For k = 2, the vertex is 

ri(- r r)r 3 (--f)r 2 (ui)r 3 (u2)r 1 (u 3 )r 2 (u4 : )r 1 (~-f)r 3 (-^) , (3.2.1) 

where r%(— 7) is the symmetrizer S2. Substituting ( |3.1.5| ) into the above, we obtain a sum of products of 
Temperley-Lieb generators ,i = 1,2,3. Mapping to Potts model is done by replacing e^'s with ( 2.5) or 
( 2.9 ), ie. convention A or B (see fi. (4) ) . 

With convention A, the corresponding operator in the Potts language has matrix elements 

W abcd (u) = ^a^(l + Q 1/2 a;i^) (q- 1/2 :e2+<W) (3- 2 - 2 ) 

a,/3,7,<5 

(q~ 1/2 x 3 + S 0S ) (l + Q 1/2 x 4 S yS ) S lc S Sd , 



where a, 6, . . . are Potts model variables taking values from 1 to Q, 

sin Ui 
sin (7 - Ui) 



and 



S a a — —Q 1 + 5, 



a/3 



7 



is the symmetrizer S2 in the Potts model language. Using the fact that 

J2 s »p = > 



the weight W can be simplified as 

W a bcd(fo, fl) = SacSbd + foS a bS c d + ^ Q 1 ^ flS a aSbaS ac S a d (3.2.3) 

a 

with 

/o = X lXM , 

h = %1 + %i + x l x 4 (Q ' + x 2 + x 3 ) . 

The three terms in the above expression have origin in the geometrical description of the T 2 model [flB]. 
They correspond to the weights of the strand configurations shown in fig.(Bl)( see appendix B ). More 
precisely, the algebra generated by the corresponding operators is identical to that obtained from the strand 
configurations with algebraic multiplication given by appending one picture to another. Further, they are 
related to the projectors Pi for i = 0, 1, 2 as 

(l)afccrf = S ac Sbd 
(Q — l)(Po)abcd — S a bS c d (3.2.5) 

and (Q — 2) (Pi) a b c d = QJ2 a SaaSb a S ac S a d — S a bS c d ■ 

Therefore we can write 

WabcdifoJl) = Wabcd + (h + Q^' 2 fl)(Q - l)(fl))«6cd + Q^^hiQ - 2)(Pl)abcd (3.2.6) 

The weight can also be rewritten in the more physical fashion 

Wabcdifo, fi) = exp [£ (Q - 1)P + K X (Q- 2)P X ] (3.2.7) 

with 



(3.2. 



e*p(fco(Q-l))-l = (h + Q- 1,2 h){Q-l) , 
exp(/d(Q-2))-l = Q- 1 /2/ 1 (Q_2) . 

This shows that the Potts model we have built is physical for 

(/o + Q- 1/2 /i)(Q-i) > -i 

and Q- 1 /2/ 1 (Q_ 2 ) > -1, 

so in particular for /o,/i > and Q > 2. It presents a mixture of ferromagnetic and antiferromagnetic 
interactions that will lead to multicritical behaviors (see section 5). It is also interesting to express W in 
terms of Kronecker delta only. It reads 

w ft n 1 - 3Q^ 1/2 /i + /o Q~ 1/2 /i - fa , x x > Q~ 1/2 fi-l fx x , rq n ^ 

W abcd (/ , Jl) = 7^5 1 ~ (Oab + Ocd)^ J\ (°ac + <>bd) (3.2.9) 

Q z Q Q 

H q (Sad + Sbc) - Q~ 1/2 fl (Sabc + 5 b cd + Sabd + S a cd) 

+ S ac Sbd + foSabScd abed 
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where 



1 if a± = a% = • • • = a ,, 
otherwise . 



The procedure of writing ( 3.2.1 ) in terms of Potts model variables is depicted in fig. (4a) which shows that 
the local interaction W abc d(fa, fi) i n Q3.2.S ) corresponds to the fundamental block Q2 which is a square with 
abed located at the four corners. Note that the interaction between the four sites are rather complicated 



as given in ( 3.2.9 ). 

On the other hand, working with convention B, we have the Potts model weight 



bcdV 



l + Q 1/2 x 3 S aa ) (q- 1/2 



(3.2.10) 



X4 



where 



S'ab 



1 



is the g-symmetrizer, which in this representation simply constrains neighboring sites to have different colors, 
in which case, it has value one. The weight can likewise be written in terms of projectors as ( 3.2.6) but 
with the projectors now given by 



(l)abc(i — S a dS 
(Q - l)(P )abcd = S' 

{Q-2){Px) ahcd = (l-S bc )S' 



(3.2.11) 



where 



<5 — S abS acS bdS cd 

It is then instructive to factor out S in the weight leaving 

W abcd (fo, h) = S'(Q- 1/2 h + S ad + f S bc ) 



(3.2.12) 



this is in fact the Potts model considered in [Q. In this form the roles of /1 and /o are more transparent; fo 
can be perceived as the parameter which controls the anisotropy while f\ controls the four sites interaction 
induced by the constraint S . This constraint imposed by the symmetrizers on the sites is however nontrivial 
as can be seen in the expansion 



W abcd (fo,h) = Q- 1/2 fi-Q- 1/2 h(Sab + 8 cd + S ac + 5 bd ) + 5 ad + f S bc 

+ (V 1/2 /l - fo) {Sabc + Sbcd) + (V 1/2 /l - l) (*«W + Sacd) 
+ Q- 1/2 fl {SacSbd + $abScd) + (fo - ZQ~ 1/2 h + l) S ab 



(3.2.13) 



ibed 



The mapping to Potts model using convention B is also shown in fig. (4b). The weight (3.2.10 ) corresponds in 
the figure to a collection of vertical and horizontal links. Notice that the lattice sites on the top (respectively 
bottom) rows are identified since they have no Ti vertex between them, and therefore carry the same color 
a ( or d ). This gives rise to the 45° rotated square on the rhs of fig (4b.). This rotated square with the 
variable a summed over is the fundamental block Q 2 and corresponds to the weight ( 3.2.13). 
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convention A 



4a 
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mapping to 
Potts model 



summing 

internal 

variables 



fused vertex 



convention B 




4b 



decomposition 



fundamental blocks 



Figure(4) The construction of the fundamental blocks 1*2 vertex is first decomposed into Ti vertices 

which occupy raw L= X, 2 and 3, the variables —7, ui, ■ ■ ■ shown next to the vertices are parameters 



that occur in ( 3.2.1 



). When convention A ( B ) is used, vertices on row i = 1, 3 are mapped to 
horizontal ( vertical ) links, and that on row i = 2 are mapped to vertical ( horizontal ) links. These 
two sets of links give rise to the fundamental blocks upon summing the internal sites variables 
labelled by a, /3, ■ ■ ■. 

The above discussion shows that the Potts model repres entati on for the T2 vertex model i s achie ved by 
replacing the spin-1 vertex either by G2 that corresponds to ( 3.2.9| ) or Q 2 that corresponds to ( [3.2.13 ). This 
is an extension of the spin-1/2 case where the Ti vertex is replaced by a horizontal Q\ or vertical Q 1 link. 

For k = 3 the vertex is 

S3S 3 r3(ui)r 4 (u2)r 5 {u 3 )r 2 {u4)r 3 (u 5 )r4 : (u e )ri(uj)r2{us)r 3 {u 9 )S 3 S 3 (3.2.14) 
with the symmetrizer S 3 given by 

S 3 = n.(-7)r2 (-27)7-1 (-7) 
S 3 = r 5 (-<y)r 4 (-2^)r 5 (-j) 

The corresponding Potts model Boltzmann weight obtained by convention A has expression 

W(u)aebcfd = ^ S aea Sb e p (S aJ 5 e fSpS + Q^ 1 ^ 2 foS a pS~fS + /i<5q^7<5 + hSa^S ps^j Sjf c Ssfd (3.2.16) 

where 



(3.2.15) 



S a bc — ~ 



1 



(1 - S ab - S bc - (Q - l)S ac + QSabc) 



Q-l 

is the contribution from the symmetrizer and 

fo = nLi^ 

/1 = xix 2 X4[x 5 + x 9 + £52:9 (Q 1 / 2 + x 6 + x s ) + xsXgiQ 1 / 2 + x 7 ) + x 6 x 9 (Q 1/2 + x 3 )] 

+ [Xl + X 5 + XiX 5 (Q 1 / 2 + X 2 + X4) + XiXiiQ 1 / 2 + x-j) + xix 2 (Q 1 ^ 2 + x 3 )]x 6 x$xg 
+ XlX2X4[(Q 1 ^ 2 + X5)(Q 1 ^ 2 + X 3 + X7) + X 3 X7}xeX S X9 + XlXg(x2X8 + X4X&) 

h = x i + x 9 + xiXglQ 1 / 2 + £ 4 + x s + XiXsiQ 1 ^ 2 + x 7 ) + x 2 + x 6 + x 2 x e (Q 1 / 2 + x 3 )} 
+ x 5 [l + xi(g 1/2 + x 2 + x 4 )][l + xgiQ 1 / 2 +x 6 + x s )] 



(3.2.17) 
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with Xi's defined similarly as those in the k = 2 case. The four terms in (3.2.16) are again associated to 
the strand configurations in the geometrical interpretation of the T3 vertex model as shown in fig.(B2). The 
weight has complicated dependence on a, 6, •••,/, it involves all possible interactions among the six sites. 
The explicit expression is given in appendix A. The construction is illustrated in fig. (5a). The top two figures 
on the rhs show that when the variables labelled by Greek letters are summed over, the resulting figure is 
C/3, a hexagonal plaquette with a ■ ■ ■ d occupying the six corners. It corresponds to the weight W ae b c fd which 
depends only on the variables a ■ ■ ■ d. 

Alternatively, one could map the T 3 vertex to the Potts model using convention B. This gives the 
Boltzmann weight 



W (u)aebcfd 



E 



Sfc 



(3.2.18) 



where fa, f± and fa are similarly defined in terms of ut as in ( 3.2.17 ), and they are decorated with tilde 
for reason that will be clear in the next few subsections. Graphically, the Potts model weig 
corresponds to a hexagon Q 3 which differ from Q 3 by 90° rotation as shown in the rhs of fig. (5b). Note that 
in this case, we have also considered the five sites e ( / ) on the top ( bottom ) row as a single site for the 
same reason as in the G 2 case. 
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L ^ J 
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fundamental blocks 



5b 



Figure(5) The construction of the fundamental blocks . The procedure is identical to the Y2 case, the Ti 
vertices occupy row i = 1, — , 5, and contain more parameters —7, u\, ■ ■ • ,ug. The two conventions 
give rise to the fundamental blocks having the same geometrical shape but differ in orientation. 

Hitherto, we have demonstrated the construction of the Potts model for k = 1 ( section 2 ), 2 and 3. 
Unlike the k = 1 case where the Ti vertex is mapped to a Potts model link (horizontal or vertical), the 
T2 and T3 vertices are mapped to a set of horizontal and vertical links. These Potts models still reside 

For 



on square lattices and have inhomogeneous nearest neighbor interaction such as that given in ( 3.2.2 ) 
the Boltzmann weight, however, it is more natural to sum over the variables associated with the internal 
sites, such as that labelled by Greek letters in ( 3.2.2 ) or fig. (5a) and work with less number of variables 
which are associated with the sites on the boundary of the fundamental block. In this case, the interactions 
are no longer restricted to nearest neighbors as can be seen, for example, in ( |3.2.S ). The corresponding 
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fundamental blocks which depend on less Potts variables are more appropriately regarded as the plaqucttes 
shown in the rhs of figs. (4a), (4b), (5a) and (5b). 

These three members of the family of Potts models appear to be quite different from one another, 
nevertheless we will show they have common symmetry properties . For arbitrary k the Boltzmann weight 
can in principle be written down following the procedure outlined for k = 2 and 3. It is expected to be quite 
complicated. It depends on k parameters and can be regarded as an operator that acts on k ( convention A 
) or k + 1 ( convention B ) Potts variables. In what follows, we shall mainly discuss the geometrical shape 
of the fundamental block Qy. ( Q k ) for arbitrary k. 

For k odd, the Tk vertex maps either to a hexagonal plaquette Qk with four slanted edges and two 
horizontal edges by convention A or to a hexagonal plaquette Q k that has four slanted edges and two vertical 
edges by convention B. The number of lattice sites on each slanted edge is equal to (k + l)/2, while there are 
two sites on the horizontal and vertical edges (see fig. (6)). The case of k — 1 is a degenerate situation where 
the slanted edge shrinks to a lattice site and the hexagon is flattened to become a vertical or horizontal link. 

It is not difficult to arrive at the shape of the fundamental block. In the Potts model language, the symmctrizcr S k , which 
consists of the vertices n, ■ ■ ■ , r i+ k-i involves (k + l)/2 sites since k is odd. These Potts sites are those residing on the slanted 
edge. The top and bottom rows of Ti vertices have odd subscripts. When convention A is used, they are mapped to horizontal 
links which becomes the two horizontal edges of the hexagon. There are 2k Potts sites on each of these horizontal edge, but 
2k — 2 of them are internal. As an example for k = 3 ( see fig. (5a) ), the four Greek letters label the internal sites. When the 
variables associated with these internal sites are summed over in the Boltzmann weight, the horizontal edge effectively carries 
two sites. In the k = 3 case, they are labelled by a c and b d. When convention B is used, the top and bottom rows of Ti 
vertices are mapped to vertical links. The lattice sites on the top and bottom rows are respectively identified due to the fact 
that there are no Ti vertices between them ( see for example fig. (5b) ). The Ti vertices which are midway from the top and 
bottom rows are mapped to vertical edges, which eventually form the vertical edges of the hexagon Q k . 



(fe+i)/2 sites 




Qk 

convention A 



dual 



(fc+l)/2 sites 
/ 2 sites 

convention B 



Figure(6) The fundamental blocks obtained by convention A and B for odd k model 



For k even, the fundamental block Qk obtained using convention A is an octagon with two vertical and 
horizontal edges each carrying two Potts sites, and four slanted edges with k/2 Potts sites on each of them 
(see fig. (7)). For k = 2, the slanted edge reduces to a single Potts site and the octagon becomes a square. 
When convention B is used, the fundamental block Q k is a 45° rotated square with k/2 sites on every edge. 

In both conventions, the slanted edges are originated from the symmetrizcrs S k , which acts on (fc + 2)/2 or k/2 lattice sites 
since k is even. For convention A, vertices at the top and bottom rows are associated with vertical edges, while vertices midway 
from the top and bottom rows are associated with horizontal edges, these edges are the two vertical and horizontal ones of the 
octagon. For convention B, all Potts sites on the top and bottom rows are respectively identified for the same reason given in 
the odd k case, they therefore become the top and bottom corners of the 45° rotated square. The other two corners of the 
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rotated square come from vertices midway from the top and bottom rows, which are associated with horizontal edges. Finally 
in arriving at the shape of the fundamental blocks, we have assumed that internal sites variables are summed over. 



(fc + 2)/2 sites 

<'//> 



2 sites 



convention B 



Qk 

convention A 



fc/2 sites 



Figure(7) The fundamental blocks obtained by convention A and B for even k model 



3.3 The patching 

The previous subsection deals with the construction of the local weight ^ of the Potts model and also 
the fundamental block Q^, '. These fundamental blocks are the building block of the lattice just as in the Fi 
Potts model the lattice is constructed from vertical and horizontal links, which, in our notation, are Q\ and 
Q 1 . Since there are two conventions ( A or B ) to be used in getting the Potts model, two lattices can be 
constructed. They are denoted as C and C . Detailed analysis shows a splitting between k even and k odd 
cass. 

For k even, we use the T2 vertex model as an example. Beginning with the vertex model lattice, we 
replace each of the T2 vertices by the fundamental block. For convention A, this gives rise to a Potts model 
lattice C which resembles a square check board [[L9|. First, each T2 vertex is replaced by the fused vertex 



( 3.2.1 ) shown in fig. (2), and subsequently this fused vertex is mapped using convention A to the square Qi 
which has only lattice sites on the four corners. Neighboring C/2's are connected along the NE-SW or NW-SE 
direction by sharing a lattice site on their common corner. For convention B, the lattice £ is a 45° rotated 
square lattice. This is obtained by replacing each fused vertex with the 45° rotated square plaquette Q 2 and 
glueing neighboring Q 2 along the slanted edge. In this case, C and C are both square lattices, however, the 
former is check board like and only alternate squares are given Boltzmann weight W, while the later has all 
squares associated with the weight W . See fig. (8) for the construction of the two lattices. 
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Figure(8) Construction of the T2 Potts model lattices where the procedure is illustrated using five neighboring T2 vertices. 

It is straight forward to extend the above to arbitrary even k. For convention A, we replace each of 
the Tk vertices by the octagon Qj~, and neighboring octagons are connected along the NE-SW or NW-SE 
direction by glueing together pairs of slanted edges. The resulting lattice C is shown in fig. (9a) ( for k = 4 ) 
where the filled circles denote the lattice sites. For convention B, the Tk vertices are mapped to 45° rotated 
square plaquettes Q k , the resulting lattice £ is a 45° rotated square lattice just like the k = 2 lattice C , 
however, there are {k + 2)/2 lattice sites on every edge (see fig. (9b) for k = 4). 




c k=4 c 

convention A convention B 

Figure(9a) The r 4 Potts model Figure(9b) The r 4 Potts model 
lattice from convention A lattice from convention B 

For k odd, we use the T3 model as an example. Recall that £3 and C/ 3 differ by 90° rotation, we expect 
therefore the lattices L and L to have similar symmetry properties. 
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convention A 




10b 



Figure(lO) The two possible patching of the fundamental blocks from the ncigboring 
vertices Vi and V2. The top figure on the rhs uses convention A for the 
neighboring vertices, the resulting Potts lattices are not compatible. 
The bottom figure shows the correct mapping where different conventions 
are used on neighboring vertices. 



As before we first replace each Fs vertex by the fused vertex ( 3.2.1 ) shown in fig. (2) and consider two neighboring fused 
vertices connected along the NW-SE direction. Let us denote them as Vi and V2 respectively (see fig. (10a)) where Vi is on the 
upper left corner of V2 ■ The Ti vertices that belong to Vi are n, • • • , rg and that belong to V2 are t±, ■ ■ ■ , rg. The symmetrizers 
Ss's that connect Vi and V2 occupy rows i = 4 and 5, ie. they are made out of r± and r$. Suppose Vi is mapped to the 
Potts model fundamental block Q3 using convention A, the Ti vertices r± and r§ will be mapped respectively to a vertical and 
a horizontal link. However, this implies that the Ti vertices on the top row of V2, which are T4, have to be mapped accordingly 
to vertical links so that the two set of Potts model links obtained from Vi and V2 are compatible. In other words, convention 
B has to be used on V2 which replaces V2 by G 3 . This point is illustrated in fig. (10b). 

Hence, neighboring T3 vertices are to be mapped to Potts model links using different conventions. The 
resulting lattice C has the geometry shown in the lhs of fig. (11). If the mapping of Vi is done with convention 
B, the lattice C obtained (see rhs of fig. (11)) is related to C by replacing all £3 by G3 and vice versa. This 
feature is present in all odd k models. For the I\ Potts model, vertical Q x and horizontal Q\ links are 
associated respectively to neighboring vertices connected along the NE-SW or NW-SE direction. For higher 
odd k, the Potts model lattice £ or £ have to be constructed with both Gk and Q k , for the same reason as 

in the k — 3 case, ie. the symmetrizer Sk contains even number of rows of Ti vertices. Since Q{ "* has the 
same hexagonal shape for k > 3, the lattice C and C are identical to that of the k = 3 model except that 
there are (k + l)/2 lattices on the slanted edge. 
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c c 

Figure(ll) The T3 Potts model lattices 



Graphically, it is easy to see that four neighboring blocks Q\ of the Tk Potts model can be combined in 
a natural way to form the fundamental block Q 2 k or Q 2k of the T 2 k Potts model (see fig. (12)). Nonetheless, 
this does not imply that a Tfe Potts model can also be considered, after summation over the appropriate 
variables, as T 2 k Potts model. Although the geometry is the same, the full definition of the Potts models 
implies that they are build from a vertex model with spin k/2 representation of U g su(2), and this constrains 
the form of the Boltzmann weights. Figure (12) looks like one step in real space renormalization group. 
The fact that starting from a Potts model one does not get a T 2 k suggests that the Tfc Potts models for 
different values of k belong a priori to different universality classes. 




Figure(12) Relation between (Jfc and <j2k- The cases of k = 2, 3 are illustrated 
which shows that four neighboring have the geometry of a §2k- 



3.4 Isotropy of the Potts model 

We consider now some symmetries of these Potts models. First, it is clear that for arbitrary k the weight 
W is invariant under reflection of the the fundamental block about a horizontal or vertical line drawn across 
its center. As an example, the weight W of the k = 2 model is invariant under the following reflections; 



reflection about vertical line 
reflection about horizontal line 



(abed) — 


-> (cdab) 


g 2 


(abed) — 


-> (acbd) 




(abed) — 


-> (bade) 




(abed) — 


-* (dbca) 





(3.4.1) 



1G 



where a, - ■ ■ ,d are Potts model variables associated to the lattice sites on the boundary of Q2 and Q 2 as given 
in figs. (4a) and (4b). 

In addition, for k even, the fundamental block has 90° rotational symmetry about its center. The weight 
associated to the fundamental block however does not remain unchanged under such a rotation. We first 
consider the k = 2 model. The 90° rotation is given, using the notation of figs. (4a) and (4b), by 

(abed) — >(cadb) (3.4.2) 

for both G2 and Q 2 - Under rotation, the corresponding weights undergo the transformation; 

W(f ,fl) — /0W90°(/ -\./c77l), r „ 4 „x 

W'(fo,fi) — /oWU/o-S/o" 1 /!), { "' 

and for fo = 1 they are invariant. We shall refer to this symmetry as face isotropy. Face isotropy is 
also present in higher even k models. To examine the behavior of the weight under this rotation, we recall 
that the underlying Tk vertex has k independent parameters which can be taken as the coefficients of the 
projectors Pj ;j = 1, • • • , k. However, a more convenient choice for our purpose is provided by the weights 
of the k + 1 strand configurations from the loop model realization of the r& vertex model as explained in 
appendix B. This set of configurations for k even contains a special element which is invariant under 90° 
rotation (for the case of k = 2, it is given by the last picture in fig.(Bl)). The others group naturally into 
k/2 pairs such that configurations that belong to the same pair differ only by 90° rotation (again, if we refer 
to fig.(Bl), the first and second pictures belong to the same pair). The parameters of the Boltzmann weight 
of the Tfe Potts model can now be taken as the weights of these strand configurations. They are denoted as 
(fo, 1), (/i,/fe-i), • • • , (/(fc-2)/2, /(fe+2)/2) and f k /2 where the brackets enclose the weights of the paired up 
configurations. Rotation of Qk or Q k about its center corresponds to rotation of the strands configuration 
and is thus given by the following transformation of the parameters; 

f0 > f0 1 j-g ^ ^ 

fifo > fk-ifo ;i = 1, . . . , k — 1, 

which is obtained by rescaling the weight ^ by fo and using the fact that configurations that belong to 
the same pair are interchanged by the 90° rotation. From the above it is clear that face isotropy is achieved 
when 

fi = fii ;i=l,...,(fc-2)/2, (3A5) 

which reduces the number of parameters to k/2. 

For k odd, there is no face isotropy since the hexagonal fundamental blocks are not invariant under 90° 
rotation. 

Next, we examine the behavior of the Potts model under 90° rotation of the lattice and shall refer to 
such a symmetry as lattice isotropy. We shall study the staggered case where weights of Gk that are 
neighboring to each other have independent sets of parameters, so the model depends on 2k parameters (in 
the following, staggered and homogeneous refer implicitely to the underlying vertex model) . Lattice isotropy 
can then exist in both even and odd k models. In fact, eventhough the lattice structures are very different 
between the odd and even k models as pointed out earlier, their parameters basically behave the same under 
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rotation of lattice. As an example, we first consider the staggered T3 Potts model. Recall that the lattice 
is constructed out of Q3 and Q 3 , and since the model is staggered, the two sets of parameters /o, /1, /2 and 
/oj/ij/2 that belong respectively to the weights ( 3.2.16), ( 3.2.18) of C/3's and Q 3 's are considered to be 
independent. To examine the transformation of the weights, we first rewrite them as 



W (/o, fl, j-l)aebcfd 



Q fo E SbeaSdfpiSapS-fgSef + Q 1 ^ f 1 S a -f5 f}5 + hfo 1 SaypS (3.4.6) 
a, (3, 7,6 

+ Q 1 ^ 2 /l/o 1 3a/3S~ ( s)S~ ( eaSsfc , 
W(fo, fl, h)aebcfd = Q _1 ' 2 ' fo z2 S aea Sbep(5 a p8 1 S + f 1 Q 1 ^ S a j8 e fS pg + fif^SaySpS (3.4.7) 



E 



+ /l/o 1 Q 1 ^ 2 $at3' 1 s)S 1 fcS$fd ■ 

After the rotation of the lattice, ( 3.4.6|) and ( 3.4.7|) correspond to weights of Q 3 and Q 3 respectively, we 
therefore compare ( 3.4.6| ) with ( 3.2.16j ) and ( 3.4.7) with ( i. 2. IS ). This implies the following relation 
between the partition functions before and after the rotation 



Z{fo, flifi] fo, flif'l) — (/o/o) Ar -^90°(/o 1 1/2/0 \/l/o i '/o i '/2/o i ;/l/o ^ ! 



(3.4.8) 



where TV denotes the total number of underlying vertices of the model. Clearly, lattice isotropy is given 

by 



fo — fo 5 

f f~ l/ 2 f ? — I/ 2 

Jl/O — J2/ 



(3.4.9) 

£ i? — 1/2 ? 1/2 

/2J0 — J1J0 

It is straight forward to extend the study of the lattice symmetry to higher k models. For k odd, we 
again adopt the parametrization offered by the weights of the strand configurations. However, there 
are now two sets of identical strands configurations from neighboring Qk and Q k , the parameters are 
(/o,l),(/i,/fe-i),-.-,(/(fc-i)/2,/(fc+i)/2) and (/o,l), (/i,/fc-i),---,(/(fc-i)/2,/(fc+i)/2) where the first and 
second sets belong respectively to the weights of Qk and Q k . Parameters are paired up according to the 
same criterion as before. Notice that in this case there is no strand configuration which is invariant under 
90° rotation. For the k = 3, these are precisely the parameters that appear in ( |3.2.16| ) and( |3.2.T8| ). Since 
rotation of the lattice by 90° turns Qk into Q k and vice versa, and also interchanges configurations that 
belong to the same pair, we have the following mapping of the parameters 



fo 

f f -V2 
JiJO 



fl 1 , 

f f-1/2 
J k—ijQ 



; i = 1, 



.,fc-l 



Thus, the model with lattice isotropy has k parameters where 



fo 

fifo 



1/2 



= fo 1 

= fk—ifo 



-1/2 



(3.4.10) 



(3.4.11) 



For even k, similar situation occurs. We again use the weights of the strand configurations as parameters, 
which read (/ , 1), (fi,fk-i), ■ ■ • , (f(k-2)/2, /(fe+2)/2), fk/2 and (/ , 1), (/1, fk-i), (/(fe-2)/2, f(k+i)/i), fk/2- 
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These two sets of parameters again correspond respectively to the weights of neighboring Gk's. The only 
difference between this model and the previous one is that the entire lattice is constructed out of Qk or Q k 
exclusively. Under 90° rotation of the lattice, the parameters transform as 

t f-fy / V-V2 ■ _ i i. 1 (3.4.12) 

JijQ < * Jk-ijQ ,1 — 1, . . . , K 1 , 

for similar reasons as in the odd k case. The Potts model with lattice isotropy is given by 

f f-l/2 _ f 7-1/2 ■ _ 1 , _ -, (3.4.13) 

J'i/O — Jk-iJo ,1 — 1, . . . , K 1 , 

and has fc independent parameters. 

Conditions ( 3.4.1l| ) and ( 3.4.13 ) both define staggered models in general. 



The above transformations of the parameters show conversely that the homogeneous model where 

a : h-i >3AM) 

is in general not invariant under the rotation of the lattice except when face isotropy is present in every Qk 's 
weight and this applies only to even k models. 

3.5 Self duality 



In the construction of the Tk Potts model, there exist two possible conventions (A or B) that can be 
used. These two choices give rise to two lattices. For k odd, the lattices C and C have identical structure 
due to the fact that both Qk and Q k have to be used, while for even k, C and C arc different. Let us discuss 
the relation between the models obtained by conventions A and B in more details. We consider a general 
staggered Potts model. We first look at the Y\ Potts model to examine the origin of this two choices. We 
have, in the notation of this section, neighboring T\ vertices X2i~\ and X 2 i given by 

= 1 + /oeai-i , fQc-n 
X 2l = l + / e 2j [6 - bA > 

where /o and /o are independent (they should be identified with x^ 1 ,X2 defined in section 2). Mapping to 
Potts model using convention A, we replace A^-i and Xn respectively by horizontal and vertical links. The 
vertices in terms of Potts model variables read 

(A 2l -i) Q 6 = S ab + foQ 1 / 2 , , . 

(x 2l ) bc = i + /oQ 1/2 ^c, [ ' 

where a b and b c are the sites on the horizontal and vertical links. After duality transformation |Q X21-1 
and Xn are instead associated to a vertical and horizontal link, which corresponds to convention B 

(X2i-i)bc = ^ + foQ 1 ^ 2 8bc, (3 5 3) 
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where we stick to the convention that a b (be) are sites on the horizontal (vertical) link. This implies, 
following fusion, that for arbitrary k, the Potts models obtained by conventions A and B are related by 
duality transformation. 

Combining the above argument with result of subsection 3.3 on the structure of the lattices C and C , 
we see that the question of self duality does not arise for the family of even k Potts models since the lattice 
C and its dual C are not the same. But for k odd, self duality can occur since both C and C have the same 
struct ure. T he self duality condition can be easily identified for the T\ Potts model by comparing ( 3.5.2| ) 



with ( |3.5.3| ). We see that duality transformation amounts to 

fa «— ► fo , (3.5.4) 



which is equivalent to ( 2.11 ). For higher odd k models, similar transformation of the parameters can be 
deduced by using the fact that duality map is equivalent to interchanging convention A with convention B. 
Thus using the same set of parameters as in the previous subsection, the duality transformation is given by 

15 y ^~ ' (3 5 5) 

ft < > fi ; i = 1 , . . . , fc — 1 . 

And self duality is given by the condition 

* : jh«-i (3 - M) 

which means that the model becomes homogeneous with k parameters. If the Tk Potts model has lattice 
isotropy to begin with, then self duality map is given by 

,°.-i/2 {° , , , (3.5.7) 

fifo < — > fk-ifo i = l,...,k- 1 . 



where the relation ( [3.4.11 ) has been used. Thus imposing the conditions of self duality and lattice isotropy, 

(3.5.8) 



the number of parameters reduces to (fe — l)/2 since we now have 

fo = fa =1 



fi = fk-i = fi = fk-i ;» = l,...,(fc-l)/2. 
Recall that in the case k = 1 these conditions determined completely the Potts model interaction 

3.6 Integrability 



We conclude this section by discussing some special cases of the Potts models which are known to be 
integrable. For simplicity, we restrict ourselves to the homogeneous models. In this restricted class, the Pi 
distinguishes itself from the rest in that it is integrable for any fo. The corresponding R matrix can be 
written as 

- , s sin u 

R(u) = l + — -e, 3.6.1 

sin (7 — u) 
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where 



/o = 



sinu 



sin(7 — u) 

For higher k, only subsets of the full parameter space are integrable. A standard such case is obtained from 
the special choice in ( 3.1.4| ) 

(i - l)fc < j < ik < k 2 , 
\ i = 1, . . . , k . 



3 

Uik 



1) 7 



(3.6.2) 



where u is the spectral parameter of the R matrix . We shall refer to this as the JB integrable line as this 
is constructed by Jimbo in Q . The R matrix for this Tk model can also be written as a linear combination 
of projectors; 



R(u) = P J=k + 



q 



21, 



1 - y 2 q 



2 r ,2k 1 



q 2k y 2 



q 



2(fc-l) 



y - q 



1 — y 2 q 2k 1 — y 2 q 2 ( k 1) 1 — y 2 q 



_ „2 ' 



(3.6.3) 



where y — exp(— iu) is the multiplicative spectral parameter. The corresponding Potts model is also inte- 
grable since the fact that Yang Baxter equation is satisfied is expressed algebraically, without reference to a 
particular representation. In terms of the parameters /i's of the T2 and 1^ Potts models the JB integrable 
line is given by the following formula 



sin(u) sin(7 + u) 



sin(27 — u) sin(7 — u) ' 



_ Q 



1/2 



sin u 



(3.6.4) 



11 sin(27-w) ' 

, sum sin(u + 7) sin(u + 27) 

sin(7 — u) sin(27 — u) sin(37 — u) ' 
, (Q — 1) sintt sin(u + 7) 

1 sin(37 — u) sin(27 — u) ' 
» _ (Q — 1) sin it 

sin(37 — u) 

In terms of the spectral parameter u, rotation of the lattice by 90° is given by 

u — > 7 — u . 

This transformation of u actually applies to all the models on the JB integrable line. The condition for 
lattice isotropy /o = 1 is therefore equivalent to u = 7/2, and the parameters /i's read 



1 



Q 



fo = 
fo = 1 

h = h = 



(3.6.5) 



(Q-i) 

Q+Q 1 ' 2 -!' 



Besides the JB integrable line there is another integrable line in all Tk models, which is related to 
the Temperley-Lieb algebra. For arbitrary k, since the projector (fc + l) q Po satisfies the Temperley-Lieb 
algebra |ft7| with 

e = (fc + l) 9 P 



and 



(k + l) q e 



(3.6.6) 
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the R matrix 

R{u) = 1 + —. (k + 1), Jb (3.6.7) 

sm(7 — u) 

satisfies the Yang-Baxter equation with u being the spectral parameter. We shall refer to this line as the 
TL integrable line. Obv iously , for k — 1, the TL and JB lines coincide. The Potts model which corresponds 
to the above R matrix ( 3.6.7j ) has parameters given by 

sin it 



t 

J0 sin(7 - u) ' (3.6.8) 
fi = ;i = 1, 1 . 

For k > 2, these two lines JB and TL do not exhaust all the integrable cases. As an example, for the T2 
model, there is an additional integrable line given by 

„ sin(w) cos(7 — u) 

JO — - ~: 



h = 



sm(27 — u) cos(37 — u) 
sinu 



1/2^, ' (3-6.9) 



sin(27 — u) ' 

which is related to the Izergin-Korepin model |l6f| . In this case the model with lattice isotropy is given by 

u = 3 7 /2 + tt/4 (3.6.10) 

or 



fo = 1, 

h = -(i + q-^ + iiq-q- 1 )) • 



4 The Spontaneous Magnetization 
4.1 Definitions 



The spontaneous magnetization of the usual (ri) Potts model on the first order transition line has been 
computed exactly by Baxter ]20]| using the method of corner transfer matrix (CTM). We will show in this 
section that using the fusion procedure, spontaneous magnetizations of the homogeneous Tfc Potts model 
can similarly be computed. As for k = 1 the model is expected to undergo a first order phase transition 
along the JB integrable line for Q > 4 p|,p4|. In what follows, we shall restrict Q to this range and replace 
7 by — iX with A being real, thus 

v /Q = 2coshA (4.1.1) 

and work exclusively on the JB integrable line. 

To define the spontaneous magnetization for the fused Potts model, we generalize the work of [^0| and 
consider some fundamental block or Q k ) sufficiently remote from the boundary and fix a site a which 
we shall refer to as central site. The spontaneous magnetization is then defined as 

M= Q(^ , 1 )-i (412) 
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where 

(^o,i)=^Po i usE^,in^' ) )- 
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Figure(13) The geometry of the Potts model lattice; a,b,c are the face variables, 
the central site <r is taken in the figure to be the filled circle at the center, 
and the spin variables are denoted by the ai where i = 1 and TV are shown. 

For the sake of computation, we consider the underlying Tk vertex model lattice to be an I x I square as 
in fig. (13) where we show the Tk( k odd in this case ) vertices and the respective fundamental blocks Gk, the 
aj's denote the spin-| arrows which can have k + 1 states §,§ — 1, . . . , — f, and i takes values from 1 to N 
with N taken to be odd always. As for k = 1 the boundary conditions are conveniently defined in terms of 
the vertex degrees of freedom: we require that the spin arrows along the perimeter all have the same state. 
This provides k + 1 boundary conditions of which we assume (more later) that they select different phases 
of the Tk Potts model. 

Recall that each spin--| arrow can be regarded from the fusion point of view as being made up of k spin-i 
arrows. These spin-i arrows form surrounding polygons for the Potts model links, in particular, due to the 
boundary condition, the perimeter can be viewed as k polygons enclosing the Potts model. We then define 
the central site a to be connected to the boundary if there is no spin-i surrounding polygon enclosing it 
other than those from the perimeter. With this definition, we can now relate the spontaneous magnetization 
to the percolation probability P defined as the probability of a being connected to the boundary. It is easy 
to see that 

(<W) =Q~ 1 (1-P) + P (4.1.3) 

and 

M = P . (4.1.4) 

As in the case of Ti Potts model, the percolation probability can be expressed in terms of variables of 
the vertex model. Referring to fig. (13), we define the following quantity pl| 

S(a) = e -(™+m<*i+-+<*N-i) ( 4L5 ) 
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and its expectation value 

(S(a)) = Z- C ] tcx £ S(a) J] (vertex weight) . (4.1.6) 
M 

Recall that closed loop formed by the spin-i surrounding polygon has orientation given by that of the spin 
arrow, and it acquires a weight (e A ) e _A when the direction is (anti-)clockwise. Writing the tVs in S(a) as 
sum of spin-i states, it is clear that when er is connected to the boundary, there are as many left pointing 
spin-i arrows among ax, ... , aw-i as there as right pointing ones, giving 

S(a) = 1 . 

On the other hand, if a is not connected to the boundary, there must be some spin-i surrounding polygons 
enclosing the central site in addition to those from the boundary. Each of these polygons includes odd 
number of spin-i states of ol\, . . . , ajv-i- The total contribution of each polygon, taking into account the 
weight from its orientation, is 

e A e -(i 7r +2A)/2 _|_ c -A e (i7r+2A)/2 _ q _ 

Thus S(a) counts the number of polygon configurations for which er is connected to the boundary. We 
therefore have 

(S(a))=P = M. (4.1.7) 

Having established the above equality we could investigate directly the spectrum of the relevant corner 
transfer matrix. There is however a faster way that uses already known results for solid on solid models. 

4.2 The local height probability and spontaneous magnetization 

The mapping between Tk vertex model and sos|2^| is standard. Height variables k's are assigned to faces 
separated by the arrow spins (see fig. (14)). The heights are given integer values consistent with 

li-li-i = 20i e {k,k-2,...,-k} (4.2.1) 

where h and U-i are heights of faces separated the spin arrow Pi. Each vertex is replaced by a square face 
with heights attached to the four corners and contributes to the partition function a weight W(li-\, k, h+i, lj) 
which is set equal to the weight of the underlying vertex. 



0t h ~ 1 ft 




Time direction 



h+i 

Figure(14) Vertex and sos correspondance 
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This sos model is a special case of the fused eight- vertex restricted solid on solid model(rsos)[^3[ The 
rsos model is integrable and is equivalent to the fusion of the eight-vertex model where heights can assume 
values from 1 to L — 1 besides the above constraint ( 4.2.l| ). In the limit p5| 



L 

I, 



oo and 

oo for all i 



(4.2.2) 



such that relative heights of neighbouring faces remains unchanged, the sos model can be recovered from 
the rsos model. More precisely, the regime III of the rsos model corresponds to the Q > 4 range of the sos 
model in the above limit | 



The Local height probability defined as|23 



P(a/b,c) 



Z(a/b, c) 



j: L a ZlZ(a/b,c) 



(4.2.3) 



has been computed exactly for the rsos model in the thermodynamics limit using the method of corner 
transfer matrix p4|. In this formula, Z(a/b.c) denotes the rsos partition function with central height given 
by a, and be are heights of the faces at the boundary of the lattice that determine the state of the arrow spin 
on the perimeter (see fig. (13)). The lattice we considered has the b face separated from the a face (the central 
face) by even number of steps (see fig. (13)). In regime III, the rsos model has ground state configurations 
such that all faces separated by even steps assume the same value while all other faces assume another fixed 
value. The sos model has also the same ground state configurations since relative heights are not affected 
by the limit. A ground state is selceted by fixing the heights be at the boundary. 



Taking the above limit ( 4.2.2), the local height probability in the the thermodynamics limit has the 
expression 

({b+c)/2-a) 2 /2k C l 

(4.2.4) 



P(a/b, 



T°° nX ((b+c)/2-a)y2k C l 



where 



in 

x 



(b 
e 



2A 



b + k)/2 
c)/2-a- 



k/2 mod 2k , 



and C l m {x) is the SU(2) level-fc string function|27|], which depends on A and has the following properties 



C l 



(I 



if I m mod 2 , 



^(tt — l - y m+2fc 



= cl 



(4.2.5) 



Cl 



C 



k-l 
k—m 



k-l 



Notice that there is no dependence on the spectral parameter u, and C l m in the above formula is nonvanishing 
because a and b are on the same sublattice. 

The spontaneous magnetization can now be computed by noting that the exponent in S(a) is 



2(ai 



OiN-l, 
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and making the substitution 

S(a) = e -^+™/mb-a) 

which depends solely on the central height for given boundary condition. The expectation value therefore 
becomes 

oo 

(S(a)) = e- {X+ ™ ,2)(b ~ a) P{a/b, c) . (4.2.6) 

a=0 

Replacing in the above formula m by m + 2nk for 

m = — k + 1, — k + 2, . . . , k 

and 

b — a = n £ Z , 

the spontaneous magnetization becomes 

Ey- 2kn 2 +2mn+m 2 /2k+k/8-l/2 ( \{m-l) / 2+nk fil 
w= -, mi = m ^ I zm (4 2 7) 

WW/ V V r r 2kn 2 + (2rn-k)n+(m-k/2) 2 /2kr , l ' V 7 

where 

/ = fe,fc- 1,...,0 

specifies the various boundary conditions. The formula can be expressed as finite sums of products of elliptic 
theta functions and string functions. Putting k = 1 into the above formula, we recover Baxter's[^0| results 
for the spontaneous magnetization of the standard Potts model. For given k, we find k + 1 spontaneous 
magnetizations assoc iated to the various boundary conditions of the vertex model. However, it is not difficult 
to show using ( 4.2.5 ) that whenever I is odd, M — 0. The latter condition is equivalent to 

J k, k — 4, . . . , — k + 2 for odd k , , 

C ~\fc-2,fc-6, ...,-& + 2 for even k , ^ > 

Depending on the parity of k, there are thus (k + l)/2 or k/2 vertex boundary conditions that actually give 
rise to vanishing magnetization. M ^ for even I, ie. 

I _ ( k — 1, k — 3, . . . , for odd k (4 2 9) 

I k, k — 2, . . . , for even k 



One then has 



with 



M EL-hi ^(2mzA,^)x- 2 /W-'/ 2 (-f)(-')/ 2 ^ 

EL-Hl 3((2m - k)i\, X 2k )x (m-k/ 2 )y2k C l m 



4 for odd k 
3 for even /c 
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where the elliptic theta functions 9 U are defined as 

e 3 (u, q ) = £y a e 2n<u 

n£Z 

d 4 (u, q) = ^(-l)"g™ 2 e 2nm . 

The above expressions are distinct for the various values of /. For example when k = 2, the results for I = 0, 2 
have different behaviors in the large Q limit as shown below. With the help of the following approximation 
for the string functions pq| 



x- m2 / 2k - ml / k C\(l + 0{x 2 )) m<0 
x -m 2 /2k-m(k-i)/k C i^ + Q(a: 2 )) m > 



C\ +m { (4.2.11) 



valid for < I < k, one finds that in the limit where Q approaches infinity, the spontaneous magnetizations 
for the various boundary conditions given by even / approach unity as follows 

M X I 1 ~ 2x ~ x 2 + (x - x 2 + s k,i% 2 )Si,o for odd k , 2 ^\ 

~ 1 — 2x — x 2 + (x — x 2 )5i t o + x&i,k for even fc . ^ ' 

Another limit to consider is Q = 4 + or A = + which divides the regions of first and second order transitions 
along the integrable linep9|]. The expansion around Q — 4 can be obtained by employing the modular 
transformation formula of the string functions, 



C ™ (r) -\Afc(fc + 2)^ £ 

' l=0m =-k+l 



e i7rmm /k sin 



n(l + !)(/' + 1) 
k + 2 



ClA-l/r) (4.2.13) 



and that for the elliptic theta functions, which is standard. The parameter of the modular group in this case 
is t — HXj-K. Making the transformation and taking the limit leads to 



e -(2fe+l) 7 r 2 (Q-4)- 1 '' 2 /8(fe+2) 



2sinbT(Z + l)/2(fc + 2)j for odd k 



M ~ S e -k^(Q-i)-^/A(k+2) ( 4 ' 2 - 14 i 
sin[7r(/ + ±)/{k + 2)\ 

Hence the magnetization vanishes with a n esse ntial singularity as Q — > 4 + . A numerical calculation of 
the spontaneous magnetization given in ( 4.2.1C ) shows that for given k and Q in the domain [4, oo], the 
spontaneous magnetizations are bounded from below and above by and 1 as physically expected. Moreover 
they are ordered as follows 

k odd M > M k -i > M 2 > M fc _3 > • • • > M (fc _i)/2 , „ , ,s 

fc even M fc > Mq > M fc _ 2 > M 2 > • • • > M fe/2 . 

As in the construction of the models, we again observe a natural splitting between k odd and fc even. 
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4.3 Conjectured phase diagram 



We have computed spontaneous magnetizations in a rather formal fashion, and their meaning for k > 1 
is not completely clear. The simplest possibility is to assume that, as in the k = 1 case, different vertex 
boundary conditions correspond indeed to different phases of the Potts model. We shall then refer to the 
cases where M = as disordered phases, although such phases may well have for instance antiferromagnetic 
order. Similarly ordered phase refers to M / 0. The order (4.2.15) suggests that such phases can be 
characterized by their degree of spin alignment. 

We can then make some conjectures about the structure of the phase diagram of the family of Tf. Potts 
model in the neighborhood of the JB integrable line ( 3.6. 2| ) for Q > 4. We shall consider only the staggered 
Tk Potts model with lattice isotropy. This Potts model has k parameters ff, i — 0, • • • , k — 1 which are 
introduced in the previous section. Recall that the JB integrable model is homogeneous and satisfies the 



condition ( 3.6.2 ), therefore the requirment of lattice isotropy fixes the spectral parameter u to be 7/2 and 
the JB integrable model corresponds to a point (denoted as Pjb in the sequel) in the k dimensional parameter 
space. We now wish to build, for fixed Q > 4, the phase diagram in the neighborhood of this integrable 
point, where the above calculation of spontaneous magnetization is performed. We expect that there are 
k + 1 distinct phases, while k parameters are at our disposal. As for k — 1 we expect that these k + 1 phases 
coexist only at Pjb- This implies that the phase diagram around Pjb has the topology of the dual of a 
fc-simplex where each phase has a common boundary with any other phases. The boundaries that separate 
the phases can be deduced from the symmetry properties of the Potts model. The discussion again split into 
two cases; k odd and k even. For simplicity we just discuss two examples. 

We start by the T3 Potts model, which has two ordered and two disordered phases. The model has a 
duality transformation given by 



fo 

and fo 1/2 fi 



fo ' 
fo h 



(4.3.1) 
(4.3.2) 



We expect that duality still interchanges respectively ordered and disordered phases, and therefore that the 
boundaries of these four phases are invariant surfaces of the duality map. These surfaces have to be given 

by 



and 



Bo 
Bi 

B12 



/o-l 
/1-/2 

F(fa 1/2 fiJo 1/2 h) 



(4.3.3) 



1/2 1/2 

where the unknown function F depends only on f Q fi and f Q J2, and satisfies the following conditions 



and 



F(fo 1/2 fiJo 1/2 h) 
F(fo 1/2 hJo 1/2 h) 



= F(f- 1/2 f 2 J 
= at P JB 



-1/2 



h) 
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fo 1/2 P. 



Figure(15) Phase diagram of the r% Potts model in the neighborhood of the 

JB integrable line point. The three surfaces fo = 1, /i = f'l and F = 
divide the four phases o, d, o' ,d' . The phase diagram has the topology 
that the four phases meet only at the integrable point Pjb- 

Incorporating the fact that the phase diagram has the topology of the dual of a (degenerate) 3-simplex, we 
arrive at the phase diagram shown in fig. (15). The phases are grouped into two pairs (o,d) and (o',d') where 
o(') and d(') denote respectively the ordered and disordered phases that are exchanged under duality. Since 
spaces below and above Bo, and spaces on the left and right of Si are interchanged by duality transformation, 
Bo and B\ are therefore the boundaries that separate o from d and o' from d' respectively. The boundary 

— 1/2 

that divides these two pairs of phases is provided by the surface B\2, which is symmetric in / /i and 

— 1/2 

/ j2, and contains the integrable point P JB . 

Consider now the example of the T 2 Potts model. The parameters are fo and fi . Notice first that the 
model with face isotropy has only one parameter, /i, and thus cannot be expected to exhibit three phases in 
the neighborhood of the JB point. This means that face isotropy has to be spontaneously broken. There is no 
self duality, but simultaneous rotation of each of the fundamental blocks by 90° about its center interchanges 
ordered and disordered phases. Rotation of the fundamental blocks amounts to 



fo 



■ 



There are two lines which are invariant under this transformation; 

Bo ■ /o-l = 0, 
Bi : G(/ ,/i) - 



(4.3.4) 



where the function G satisfies 



and 



G(/o,/i) 
G(/o,/i) 

dG 

Wo 



G(fo\h) 
at Pj B 

. 



29 




Figure(16) Phase diagram the T2 Potts model in the neighborhood of the 
JB integrable point. The three phases o, o', d are separated by Bo 
and B\ , which are invariant curves under rotation of all the faces. 

The typical phase diagram is shown in fig. (16). 

Finally we want to remark that the above scenario relies upon the assumption that there is no other 
phase transition lines that bifurcate from the JB integrable line at some Q > 4. Recall that such phenomena 
occurs in the Ashkin- Teller model where the self dual line (where the underlying staggered vertex model 
becomes homogeneous) bifurcates into two phase transition lines. 

5 The T 2 Model 

While the models for Q > 4 are generally expected to be non critical, the phase diagrams for Q < 4 should 
exhibit several kinds of second order phase transitions. We shall here discuss in some details the case of 
IV We restrict to the non staggered case (which would be the two lines x 2 = 1 in the Ti case) and to the 
geometry of cylinder, ie. periodicity in time direction, with free boundary condition on the top and bottom 
rows. This ensures quantum group symmetry, which will turn out to be quite a useful ingredient. To start, 
we discuss the related one dimensional quantum spin chain. 

5.1 The quantum spin chain 

The hamiltonian can always be written as a general linear combination of the projectors as follows; 

21-1 

H=^2 ( sinw - cosw)(Q - l)P (i, i + 1) - cosw(Q - 2)P 1 (i, i + 1) (5.1.1) 

i=l 

where Pj(i,i + l)'s are projectors that project onto the irreducible spin-j representation from the tensor 
product of the spin-1 states at sites i and i + 1, and the spin chain has free boundary conditions. We have 
chosen the coefficients of the projectors for later convenience. The parameter ui takes values in [0, 2w] and q 
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is restricted to the case \q\ = 1, we define as before q = e 17 ,7 € R and introduce the parameter (5 = 71/7. 
Owing to the existence of the unitary transformation 



UPjiq)!/- 1 = Pji-q- 1 ) ;j = 0,1,2 (5.1.2) 



where 



and 















\ 


V2 





i) 


, s*= 1 


• 








V2 








-1 / 



17 = 1 + -S z <g> S* - -S* 2 ' ® S 12 ' - ~S+ a 5- 2 - is*" 2 ® S+ 2 
2 2 4 4 

+ (3S + S Z <g> (1 - 5 2 )5" + ^S'S* <g> (S 12 - 1)S*+ + aS+(l + 5+) ® S* 2 ^ 
+ a- 1 S-(l + S z )<E>S z S + ;a,(3eC 

( V2 \ 
5+ = >/2 , 5" = 
V / 

are the usual su(2) generators, it suffices to consider 7 in the domain [0, 7r/2]. The hamiltonian is in general 
not hermitian with 

H\{q) = H^q- 1 ) , (5.1.3) 

however, through relabeling of the spin sites, H^q" 1 ) can be made equivalent to Hi(q) and therefore the 
energy eigenvalues are real or occur in conjugate pairs. The hamiltonian can be written in terms of the more 
familiar su(2) spin operators as, 

2/-1 

H = {—2 cos a; cos 27 + cosw — sinw + coswcr.; + sinwof + sin 2 7[(sinw — coscl>)((7 2 — er 2 ) 
»=i 

„^2 i sin 27, . isin.7. , . 
- 2cosw(5? + H — -(sinw + cos w) of (5? - H ^-(smw + cosw) 

[a-^Sf - S? +1 ) + (5 2 - S^)^] + (sine, - cosc)(cos 7 - l)(a^af + a?^)} 

+ icoswsin2 7 (S 2 - S^) (5.1.4) 



where 



<Ti = Si ■ Si + \ = a i + a 



1 > 



Z QZ QZ 

a i — D i °i+l 

Written in this manner, this hamiltonian can therefore be regarded as a special case of the spin-1 XXZ chain 
where the boundary terms ensure U g su(2) symmetry. 

The hamiltonian reduces to that of the bilinear biquadratic spin chain with su(2) symmetry at 7 = (or 
q = 1). It has the simple expression 

21-1 

H = [cosuxTi + sinwCT 2 ] . (5.1.5) 
i=i 

This model has a nontrivial phase diagram. We summarize in fig. (17) and below certain features of the phase 
diagram M: 
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(3,3)- 
(3,3) - 



(3,3)- 



(3,3) - 




Figure(17) Phase diagram of the T2 spin chain. The various integrable lines have central charge given by 



IK 

FZ 
TL 

IK 

FZ 
TL 



12 
4(4-2) 



-1 C : 

-1 C : 



2 4(4-2) 
1 6 

-L TT— —1 



C = 1 

C = 1 
C = 1 



massive 
6(4-l) 2 



4 



S' («'-!) 



5 e [2,6) 
s e [6, 00] 

5 6 [2, 00] 
8 G [2, 6] 

<5 e [6, 00] 

8 G [2,6) 
8 G [2, 00] 
S e [2, 4] 



where 5 is related to 5 by — 2ecos -p- = 2 cos 2ZL + 1. Some of these values are indicated in the figure. 

The phase diagram is essentially divided into four regions ui G (f , ^ L ), (^f, ), (— : 



and ("2L 1L\ 

_ V2!4/>V4>2<"V2>4-' V 4 i 2 / 

depending on the ground state of the spin chain. For ~ < u> < the ground state is ferromagnetic. 
The model is integrable at uj — ^ and has (formally) central charge c = —00. At the boundaries, u> = ^ 
and ^ the symmetry is augmented from su(2) to su(3). More specifically, at ^, neighboring spins assume 
the representation (3, 3), the hamiltonian being proportional to 3-Po is related to the spin-i Heisenberg 
antiferromagnetic spin chain via the Temperley Lieb algebra. At ^ neighboring spins belong to (3,3) 
representation, the spin chain is the permutation model studied by Sutherland et a/pl|] and is found to 
have central charge equal to 2. In the interval ^ < uj < ^L, the spin chain is found using a semi classical 
approach to have vanishing magnetization but nonzero tensorial order parameter, the ground state therefore 
exhibits a "nematic order" . At ^ the spin chain has hamiltonian — Pq which differs from that at | by a 
sign, it is again in the representation (3, 3). The ground state is found to have massive excitation. The 
interval ^ < uj < ?, where we have identified the point uj = and 2ir, has antiferromagnetic ground state 
and contains the Takhtajan-Babujian model|3^] at uj = this point is solvable with c = |. In addition, 
the exact valence bond ground state @] can be constructed at uj = tan -1 -|, and the spin chain is shown to 
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have massive excitation. The vicinity of this point, which includes the Heisenberg antifcrromagnetic model 
at u> = 0, belongs to an antiferromagnetic fluid phase or disorder flat phase|3^] where there is long range 
antifcrromagnetic spin order and position disorder. For ? < u> < the ground state is dimerized, and at 
to = ?, where phase transition occurs, the model is integrable and has su(3) symmetry. 

5.2 The Integrable Lines 



The 1?2 model has mainly been studied along the integrable lines ( 3.6.4 ),( 3.6.E ), ( 3.6.8| ), they are given, 
in terms of parameters of the spin chain, as 



(FZ) 
(IK) 
(TL) 



tanoj 



tanw 



cos u> 



-3 



. 



The first two are related respectively to the and A| solutions to the Yang-Baxter equation! 



they 



have been studied first by Fateev-Zamolodchikov|35| and lzergin-Korepin|f6|. The one labelled by TL has 
hamiltonian proportional to the spin-0 projector which is known to satisfy the Tcmperlcy-Licb algcbrajl7|. 
We also want to point out that the FZ and TL lines are respectively the k = 2 element of the family of 
integrable models denoted as JB and TL in previous sections. Each of these equations gives rise to two 
lines in the lo — 7 phase diagram where the hamiltonians differ by an overall sign. In the limit 7 = 0, they 
reduce to the integrable points of the su(2) invariant bilinear biquadratic spin chain. In this section, we 
shall examine the phase diagram beginning with these integrable lines, they will serve as benchmarks for 
the understanding of the critical properties of the general phase diagram. In fig. (17), we summarized the 
features of the phase diagram. 



5.2.1 The TL case 

It is governed by the hamiltonian 

where 



H = € 



21-1 

E 

i=l 



(Q-l)P (i,i + l) 



(5.2.1) 



1 if w= 7r/2, 
-1 ifw = 37r/2 



The projectors (Q — 1)Pq satisfies the Temperley Lieb algebra ( 2.6) with 

ei = (Q- 1)P 

and 

e? = (Q - l)ei . 

The corresponding vertex model has transfer matrix that satisfies the Yang-Baxter equation, the model is 
therefore integrable. In principle, energy eigenvalues and hence the critical properties can be deduced from 
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the Bcthe anatz solution. However, we shall instead employ all that is known about the U 9 su(2) invariant 
spin-i chain||, which also has hamiltonian given by sum of Temperley Lieb generators, to understand this 
integrable case. Since the two hamiltonians are related to the same algebra, they share the same set of 
eigenvalues. On the other hand, as the representations are different, we do not expect the degeneracy to be 
identical. Also the same eigenvalue may appear in different spin sectors in the two models. To overcome 
these difficulties, we compare numerically their eigenvalues. It is worth pointing out that the hamiltonian 
has a hidden U g sl(3) symmetry [|36| where neighboring spin sites are in the (3,3) or (3,3) representation of 
the quantum group, and Po can be regarded as the operator that projects the above representation onto the 
trivial representation. 

The spin-i U g su(2) invariant spin chain with free boundary condition has hamiltonian 

21-1 

H=-J2VQ r Po{i,i + l) (5-2-2) 

i=l 

it is the extreme anisotropic limit of the self dual six vertex model ( 2.16] ), and the projector y^Q 7 Po satisfies 
the Temperley Lieb algebra with 

&i = \/~Q' Poii, i + 1) • 

The spin chain is critical for \[Q £ [0,2] (or S G [2,oo]), and the central charge depends on <5 aspsj 

The ground state energy of spin-j sector ej scales asQ 

K4 - e) ) i 

where 



^ (5.2.4) 



h . = -±1 tl (5.2.5) 



j[j(S 

5' 

and the sound velocity 

5 IT 

c = — sm —7 

^ 2 5' 

is obtained from the Bethe anatz solution ]4l[|. For q a root of unity, the central charge belongs to the minimal 
series Q with 

c = l . 6 ~r (5.2.6) 

m(m + 1) 

and m = 6 — 1. The conformal weight is given by 

\{m + l)r — ms} 2 ' 



(5.2.7) 



4m(m + 1) 
and therefore 

hj = h 1A+2j . (5.2.8) 
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For Q > 4, the spin chain is noncritical and has a massive excitation]]]]. 

The spin-i chain also has a ferromagnetic counter part[[l], whose has hamiltonian differs from (5.2.2 ) by 
an overall sign. The negated hamiltonian can in fact be obtained from ( 5.2.2 ) by rewriting the coefficient || 



-2cos(tt(1 - — )) 



and extending the domain of S' to [1,2] so that (1 - jr)^ 1 £ [2, oo] or -y/C? £ [0,2]. At y/Q 7 = 0, it 
can be shown that the eigenvalues are symmetric about zero, the hamiltonian is therefore equivalent to 
its ferromagnetic counter part. We can therefore regard the ferromagnetic counter part as an extension of 
( 5.2.2) ) where the domain of y^Q 7 is enlarged to include [-2,0] or 5 £ [1,2] as well. More importantly, it is 
f ound that for S £ [1, 2] the central charge and conformal weights are correc tly given by ( fr.2.3| ) and ( |5.2.4| ), 
( 5.2.5 ). In other words, the above results apply to the hamiltonian ( |5.2.2 ) with 6 £ [1, oo]. 

At \J~Q = —2, the spin chain is noncritical with c = — oo, and for y/Q r < —2, the ground state has 
massive excitation. 



Comparing the hamiltonians ( 5.2.1 ) and ( 5.2.2| ), we see that the two spin chain have the same set of 
eigenvalues when 

\/Q 7 = e(l-Q) ) (5.2.9) 



since the Temper ley Lieb algebra realized by them have the same y/"Q" parameter. Numerical studies of 
the eigenvalues for the spin-1 chain of finite size ( 2/ < 10 ) reveals that the spin-1 energy spectrum contains 
many crossing of eigenvalues due to the additional U 9 sl(3) symmetry In particular we find, 



-2j-l 



(5.2.10) 



always hold. By comparing the ground state energies for various spin sectors of the spin-i and spin-1 chain, 
we find, when Q and Q are related by ( |5.2.9| ), 



-1 : eh = e'l ;j>0 for Qe[0,4] 



-2.; 



e' 1 ;j>0 for Q £ [0, 2] 



(5.2.11) 



where eigenvalues with prime belong to the spin-i chain. For e = — 1, this identification implies that the 
spin-1 chain is critical for Q £ [0, 3], the central charge is given by ( 5.2.3] ) with 6 related to 5 by ( |5.2.9] ) or 



7T 2-/T 

2 cos -7 = 2 cos — + 1 . 


and it increases from -7 to 1 as Q varies from to 3. Moreover, the ground state of the spin-j sectors scales 

as 

Kel-elj) n^oo j\j(S' - 1) - 1] 

^7T 5' 

with the sound velocity given as before. Using ( 5.2.10] ), the scaling behavior of the odd spin sectors can 
also be deduced. Since Q £ [3,4] is mapped to the noncritical region y / Q 7 £ [2,3] of the spin-i chain, the 
spin-1 chain in this interval is therefore massive. 
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For e = 1, the identification is valid for Q G [0, 2] and the spin chain is critical with central charge given 
by ( 5.2.3 ) but in this case 5 and 5 are related by 



7T 2-7T 

2 cos -7 = — 1 — 2 cos — , 
o o 



hence the central charge varies from to -7. For Q £ [2,4] the mapping ( |5.2. ll| ) does not hold anymore, 
and it is not clear how to use the spin-i chain to deduce the critical properties of the spin-1 model. 



5.2.2 The FZ case 

The FZ integrable spin chain has hamiltonian 

21-1 



H = e^2[2(Q-l)P (i,i + l) + (Q-2)P x (i,i+l)} ; e = ±1 , 



it is the extreme anisotropic limit of the vertex model given in ( 3.6.2 ) for k = 2. 
the model is critical with central charge given byE3] 



For the e 



12 



5 E [2,oo) 



2 8(6-2) 

while the lowest eigenvalue of each spin-j sector we found numerically to scale as 



1(e) 



i 



where 



25 

7r sin 26 
26 



-'j.odd 



(5.2.13) 
T regime, 

(5.2.14) 
(5.2.15) 



denotes the sound velocity. For q a root of unity, 8 becomes rational, the central charge belongs to the 
superconformal series E^] where the conformal weight reads 



_ (p5-q(5-2)f 
p ' q ~ 85(5 - 2) 

Substituting p — 1 and q = 2j + 1 into the above, 



1 



P-91 



In 



2j + l 



m - 2 )j - 2 ) 

25 



(5.2.16) 



(5.2.17) 



we recover ( 5.2.15 ) except for the additional term |5 Ji0 dd ■ The spin chain is therefore related to the Neveu- 
Schwarz sector of the minimal superconformal series for q a root of unity. Furthermore, only the lowest 
eigenvalues of the even spin sectors are related simply to the primary states \h\^+2i >, while for the odd 
spin sectors they are related to G_i|/ii,t+2j > where G_i is the fermionic raising generator of the global 



superconformal group OSP(2|l). This extra factor accounts for the term |5 J)0 dd m ( 5.2. IE ). Another 
interesting phenomena occurs at the point 7 = ?. It has been noted that the numerical estimate of the 



36 



central charge is exactly zero and does not suffers from finite size correction. This point is related in fact to 
the N — 2 supersymmetric series |D|. 

The e = 1 regime of the FZ line has drastically different behavior from its e = —1 counterpart pgf] . It has 
been studied for toroidal boundary conditions where 

o± _ p ±i0c± oz _ qz 

°2l + l — L °1 ) °2l + l — °1 ■ 

The "effective" central charge depends also on tfi as 

c=l-g; 7€[0,J. (5.2.18) 

It is however well known that with appropriate value for 0, this formula gives the central charge for the free 
boundary spin chain Jil[ ftij . Indeed, putting = 27r — 27, we get 

6(8 - l) 2 

c = 1 - K ' ; 8 e [2,oo]. (5.2.19) 

This expression is the same as that for the spin-i U g su(2) invariant spin chain in the domain 8 6 [1, 2] as 

can be seen by replacing 8 by yrrr m ( 5.2.3] ). We have also verified numerically for the e = 1 FZ line that 
the ground state energy scales according to 

where 

7r sin 2(5 
^ ~~ 2ir-2S ' 

The above formula is in fact equal to ( [3.2.5 ) after the replacement 

r' t 

5 ' y~i • 



This regime is therefore in the same universality class as the spin-i chain ( 5.2.2 ) in the interval 8 G [1,2]. 



5.2.3 The IK case 

The IK integrable spin chain has hamiltonian given by 

21-1 



H = iTTm ^2 E ^ 4 -QXQ- Wit i + 1) + (3 - Q)(Q - 2)Pi] (5.2.21) 



where e = 1 corresponds to 

^ = tan- 1 (^i-^) G [0,tt] 
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and e = — 1 to 

id = tan - (- 



-3' 



The model coincides with the FZ chain at 7 = j and the TL chain at 7 = ^ . Exact Bethe anatz solution has 
been worked out for the model with toroidal boundary condition pE}]. It was found that the critical behaviors 
are classified according to the following regimes; 



regime I e = 1 c = 1 - $fe I7e(0,f) 

regimell e = -1 c = § - ^|^y I7e(f,f) 

f 2-|£ 0<2 7 

regime III e = -1 c = ^ 2 ^ r , , 2 r " ;7G (0,£) ; . 



(5.2.22) 



As in the previous case, these results can be used to obtain the central charge for the free boundary spin 
chain. We verify numerically that in regime II the central charge is given by <fi = 2j where the above formula 
becomes 



This expression is identical to ( 5.2.14| ) of the e = 1 FZ line. Moreover, numerical check of the energy 
eigenvalues shows that the ground state of each spin-j sector scales as ( 5.2. 15| ) but with the sound velocity 



given in this case, following J15| by 

2ir sin 27 cos 37 



(tt - 67VQ + Q(Q - 3) 2 



We therefore conclude that regime II is in the same universality class as the e = — 1 FZ line for 7 G (fi §)■ 
However beyond 7 = ? ie. regime III, the IK model has different critical behavior. Our numerical checks 
of the conformal weight proved inconclusive due to poor finite size convergence. On the other hand, it is 
known that the spin chain at 7 = is related to the permutation model studied by Sutherland et al where 
c = 2 and to TL model at 7 = ^ where c = 1 . It is therefore likely that regime III has 

c = 2-| (5.2.24) 


which is obtained from above with the substitution <fi — 27. As for regime I, numerical check suggests again 
that the critical properties are again different in the two regimes 

regime I e = 1 7 € (0, ^) 
regime i" e = 1 7 <E (f , f ) . 

In regime I , the model is found to be in the same universality class as the e = — 1 FZ line where the central 



charge has expression ( 5.2.19) which is obtained from the above by taking (f> — 2-k — 27, and the scaling 
behavior of the energy eigenvalues is given as in ( 5.2.2C| ) with sound velocity 



2ir sin 27 cos 37 



3(tt - 2 1 )^Q + Q{Q-Zy 
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for regime I , finite size convergence is poor and classification of the regime is uncertain. 

It is intriguing to find that the IK and FZ line in both e = ±1 have the same critical properties for 
7 € (§,§)• To elucidate this we performed further numerical study for models "in between". The sound 
velocities are not known then, and the difference of ground state energies ej — is more difficult to use 
to deduce the critical behavior ( such as ( 5.2.15| ) ) of the model. One can still study quantities that do 



not depend on the sound velocity and compare them with the FZ and IK integrable cases. For given 7 
in the shaded regions a and b in fig. (17), we found that such scaled quantities approach those common to 
the two integrable lines as I increases. Further, the ordering of levels with respect to j is the same as that 
of the integrable lines. This behavior suggets that the shaded regions are massless phases, with the same 
universality class as the integrable lines. Other indications come from the fact that on the integrable lines, 
the ground state energies of certain spin sectors coincide at special values of 7 such as 

e = -1 fel(f) = (5 2 25) 

and 

e = 1 fei(f) = 

fci(f) = 54(f), (5.2.26) 



which can be seen from ( 5.2.15 ) and ( 5.2.5| ) respectively. Numerical study shows that such crossings still 



hold in the shaded regions. Finally at 7 — j in region a, moving away from the integrable lines along the 
ijj direction, the central charge is found to be exactly zero without finite size correction as in the integrable 
cases. A study of the operator contents of the continuum theories on the FZ line suggests that the operator of 
the N = 1 supersymmetric series that correspond to perturbing spin chain from the FZ line in the u> direction 
and which respect the quantum group symmetry of the spin chain has conformal weight h^i = (8+2)/ (28—4) 
ie is irrelevant for 7 > ir/6 (region a). As for region b, it is likely that the operator of the minimal series that 
drives the perturbation from the integrable line as again weight h% i = 28—1, ie is irrelevant for 7 S [0, tt/2]. 

5.3 P2 projector and the g-deformed valence bond states 



Besides the integrable lines, the hamiltonian obtained by summing projectors P 2 (i,i + 1) deserves further 
examination. Recall that for q = 1, the ground state of vanishing energy could be exactly constructed using 
Valence Bond States JtJ. It turns out that the construction generalizes to arbitrary q. We shall again refer 
to the corresponding state as VBS. Notice that for q = e n , such a state needs not always be the ground 
state. We shall in fact observe other eigenenergies crossing as 7 deviates from zero. 

The hamiltonian considered lies along the line 

tanw = — — ;we[0,7r] (5.3.1) 
Q — 1 

with 

JV-l 

H = ^(Q-1)(Q-2)P 2 (M + 1), (5-3.2) 

i=i 



3!) 



where we have dropped an overall positive coefficient and a constant term. This restricts to the (q = 1) 
su(2) invariant model as a special case. Let us now extend the valence bond state method, used so far for 
the su(2) model, to the above hamiltonian. As a first step, we regard the spin-1 state ip a /3 as being formed 
by the q- symmetric product of two spin-i states defined as 

= 9aP<t>a ® <t>P + 9pa<\>(i ®<j>a \ OL, (3 = 1, 2 (5.3.3) 

where g a p is the matrix element of 

1 V2qi 

s/2 q+q" 1 

and (p a denotes the orthonormal basis of the spin- | states. The spin-1 state ipap is by construction symmetric 
in the two indices a, (3 and is related to the orthonormal basis |± >, |0 > as 



ipn = V2\+ >,-0 22 = \/2|- >,Vi2 = "021 = J , J O > ■ 

V ? + 9 

Note that in our notation, the two spin-^ states are labelled by 1,2, and the three spin-1 states are labelled 
by ±, 0. Consider two neighboring spin-1 states ip aai and i^^p, we construct a tensor product using two of 
the spin-i states (one from each spin-1 state) such that the total spin can only be or 1, this is achieved 
with the help of the tensor 

and the tensor product is defined as 



ft a p = ^2 ^aaieai/3iV>/3i/3 ■ (5.3.5) 
a, Pi 

In the q — 1 case, e reduces to the usual antisymmetric tensor with e 12 — 1- One can check that the resulting 
element indeed belongs to the U g su(2) spin 1 representation by expressing the above as 

On = 7 =L= (q 2 + q- 2 )i\l,l> , 

^22 - -y=L=(q 2 +q- 2 )?\l,-l> , (5.3.6) 

vq + q 

qCt^+q- 1 ^ = -2(q 2 + q- 2 )?\l,0> 

and fii2-fi 2 i = 2(q 2 + l + q- 2 )^\0,0> 

where {|1,± >, |1,0 >} and {|0, >} are respectively the orthonormal basis of the irreducible spin-1 and 
spin-0 representations, which are constructed out of two copies of spin-1 (four copies of spin 1/2) orthonormal 
bases from sites i and i + 1 as follows 

spin-0: |0,0> = (q 2 + 1 + g -2 ) - ^- 1 ! + - > |00 > +q\ - + >) 

|i,i> = (<z 2 + <r 2 H(< ? |o+>- ( r 1 | + o>) (5.3.7) 

spin-1 : |1,0> = {q 2 + q- 2 )-i (\ - + > +{q - q-^OO >) - \ + - >) 
|1,-1> = {q 2 + q- 2 yHq\ -0 > -9 _1 |0- >) • 



40 



These formulae show the important fact that the tensor product ( 5.3.5 ) satisfies 

P 2 (i,i + l)n af3 = . (5.3.8) 

Such a construction can be extended to the whole spin chain by tensoring neighboring spin-1 states with e 
giving _^ _^ 

Maf) = Y ^af3 1 £f3 ia2 '4' a2 (3 2 ---ef3 N _ iaN tp aN (5.3.9) 

ft 

ie[2,N] j&[l,N-l] 

the VBS state. It satisfies 

frng } = o (5.3.10) 

which can be checked by considering the action of individual P2(i, i + 1). 

For the chain with free boundary conditions, the indices a, (3 give rise to four states which, when expressed 
as linear combination of strings of |± >, |0>, have the characteristic that a nonzero state | + > (| — >) must 
be followed by a |0> or | — > (|+>). Thus |+ > and |— > appear alternately in the VBS and there can be 
any number of |0> between the |+> and |— >. The four states are distinguished by the following 

ilii The first nonzero states is |+> and the number of |+> states exceeds that of the |— > states by 1. 

0^2 The string has equal number of |+> and |— > states or all |0> states. 

f2|x Same as in the fii2 case. 

f2rj2 Same as in the On case with |+> replaced by |— >. 

From the lattice gas point of view (where |0 > is regarded as a vacancy), the VBS exhibits a perfect 
antifcrromagnetic order and positional disorder. 

Since Clg (^22) contains an extra |+ > (|— >), we have, 



S — ^11 

} N - -n 1 

L 22 — "22 ; 



s z ng = -ng, (5.3.11) 



and S z n? 2{21) = . 

As in the q — 1 case, it can be proved that the four states belongs to the spin-0 and spin-1 irreducible 
representations, namely, 

fin « |1,1>, 

ng cx |i,-i>, 

qQ^ + q-^g cx [1,0 > 
and tt? 2 - fig cx |0, > . 

The norm can also be computed. We define scalar products by treating q formally as a real parameter so 
that the conjugate of the raising operator is the lowering operator, and vice versa. As an example the state 

g|+> +g~ 1 |-> 
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has norm q 2 + q 2 instead of 2. With this convention, eigenstates with different eigenvalues continue to 
be orthogonal for complex q. But we loose positivity and definiteness in general. The computation is now 
done using graphical means which generalize the method of . The contraction (fi^ , ) is represented 
graphically as two parallel series of horizontal links and dots (see fig. (18)) 

Q af 3 a • • ••— • • • • • • • • « ■ • p 

Q^s 7 • • • • • « »« » • • i ■ • 6 

Figure(18) Graphical representation of the valence bond states 

where each pair of closely spaced dots represents the two spin- ^ states at each spin-1 site and the horizontal 
links represent the presence of the valence bond, ie. the e tensor. Each pair of dots has contraction only 
with that directly below (or above) it, which gives the contraction of the spin-1 states from il^g and Cl^p 
at the same sites. We first examine the one particle norm 

(ip-yS, ^ap) = K al3 (S al S/3S + SagSffy) (5.3.12) 

where 

K »P = dip + 9p a ■ 



The rhs of ( 5.3.12 ) can be represented graphically as (figs. (19a) (19b)) 

a f3 a 



X 



'8 

! 

Figure(19a) Figure(19b) 

Contraction of two spin-1 states 

and we shall refer to these two geometrical objects as the parallel and crossed vertical links respectively. 
It is now clear that ) gives 2 N possible graphs which are obtained by replacing each of the one 

particle contraction by fig. (19a) or fig. (19b). A typical graph for N = 4 looks like 

a 13! Pip~i Pi0 p p 

imn 

7 ft Aft ft/3 P S 
Figure(20) Typical graph for N=i 



which carries a weight 

5 *P K <xPi (eftft e| l/3l )^ftft ( e pp S ps) 

Pi 

and we have introduced the notation 

1 if $ = 2 , 



P = 



2 If 0=1 
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Our task now is to sum up the weights of the 2 N graphs. Unlike the q = 1 model where the sum can be 
performed using combinatoric arguments only, the q dependence of K a p and e complicates this approach 
and the sum has to be done with the help of recursion relations. 

We shall define the collection of horizontal or vertical links which are connected together as a circuit, 
thus each graph is made up of disconnected circuits. Notice that the circuits come in two different forms, 
which arc distinguished by the type of vertical links at the right and left most ends. In the above example, 
there are three disconnected circuits; The one in the middle has its rightmost vertical link formed by one 
of the parallel vertical links fig. (19a), while the circuit on the right has its rightmost vertical link given by 
fig. (19b) and left most link given by fig. (19a). We shall refer to them as circuits of type A and B respectively, 
note that in our definition, circuit of type B is characterized by the vertical links at its two ends, while circuit 
of type A by its rightmost end only. We also introduce the notion of length for these circuits, namely the 
length of a circuit is equal to half of the number of valence bond it covers. Therefore the B and A circuits 
in the example have length 2 and 1 respectively. Having established the notations, we are in a position to 
characterize the graphs. Any graph belong to one of the following types: 



graphs whose rightmost circuit is an odd length type A circuit, 
graphs whose rightmost circuit is an even length type A circuit, 
graphs whose rightmost circuit is type B, 

graph which does not contain any parallel vertical link given in fig. (21). 



In the above definitions, N denotes the size of the spin chain, c and o denote the parity of the rightmost 
circuit. It is not difficult to see that the above four cases exhaust all the possible types of graph and are 
mutually exclusive. Among the 2 N graphs, only one of them is of type g, it is made up two disconnected 
circuits running from one end of the graph to the other as given by the last figure in fig. (21). Since the first 
three types of graphs always come with the factor S ai 5ps, we explicitly separate the factor from the rest of 
the contribution. Graphically these four types of graphs have the following features (fig. (21)) 



e(JV) 



Odd 



f even f 



jOC--- Q( 



Figure(21) The four types of graphs 

The weight #), fiC), CW of N sites can be related to that of N + 1 sites as 

A o(N+l) _ .e(N) 2 R W 2 
A ap - A a p e pp a a p e pp > 

A e(N+l) .o(N) 2 

a P " ocp t (ip ' (5.3.13) 

U a p - 2^p 1 \ A a p 1 e p 1 p 1 ^l3 1 p + A a p 1 C p 1 p 1 *Pi0 + ^aP^p.p^^p) 

+ ^aa^aP^N,odd + ^a/3<W,even , 
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which correspond to the various ways of appending an additional site to the AT-site graphs (see fig. (22)) 



WTk. 

n f 



n 



A 



n cvc 



3(1(1 



ra^mc- mi 



even , 
or odd 



Figure(22) Graphical representation of the recursion relations ( [3. 3.1c ). 

The weight of can be calculated directly as 



{K ai ) N 1 Kp S e a pe lS ; N S even , 
{K ai ) N 5 a s8 lP ;iVeodd. 



(5.3.14) 



The set of recursion relations ( 5.3. 13| ) can be solved easily when they are iterated once to relate graphs 
whose length are of the same parity. This gives, for TV odd, 



.4' 



o(JV+2) _ 
a/3 

e(N+2) _ 



A 



B 



*a/3 
(JV+2) 
a/3 



= A 



a/3 
e(JV) 
a/3 



Eft « ( ft ] k m + A^k PlP + b%\k m ) + k a(i 



B 



(AO 



= Eft [ A aS? K *f> + A «Z'KW + B £t*M + Ea^^ft^ft/3 



a/3 



(5.3.15) 



where 



K af3 = £aa^a/3el 



00 



(5.3.16) 



The above may be written more compactly as the matrix equation 

G (N+2) +1= ( G (/v) + 1 ^ 1 + K)(l + K) 

where 

q(JV) = A o(JV) + A c(JV) + B (N) 

is a 2 x 2 matrix whose indices are labelled by a and j3, and from the definition, it includes 2^ — 1 graphs, 
the missing one being C^ s . 

For N even, similar calculation gives 



where 



q(n+2) +1 = ( G (/v) + 1 ^ 1 + K )(l + K) 



(5.3.17) 



^*0 e 00 
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Equations ( 5.3.16| ), (5.3.17) lead to the results 



GW = (1+K)[(1 + K)(1 + K)] ; 
GW - [(1+K)(l+K)]f -1 

Taking q = 1, the rhs of the above formulae reduce to 

3 W - 1 / 1 1 
1 1 



N € odd , 
N 6 even . 



(5.3.18) 



(5.3.19) 



which is precisely the result derived in |t]]. For arbitrary q, these formulae can be further simplified by noting 
that 

(1 + K)(1 + K) = I -^ T7 (l + (q2 + q-*)(q + q-i)i-p) 
(1 + K)(1 + K) = Iiq JL TF (l + ( 9 2 +9 - 2)(g + g -i )2 p' ) 

where 



(5.3.20) 



and 



satisfy the property 



(<?+<r 
1 

q + q- 
l 

q + q- 



— i — ? 

q q 

2 

r q 



q q 

q 2 q- 1 



p(')2 = p(') 



This relation, when combine with ( 5.3.14 ), gives the result of the contraction 



1 ( A N -(-l) N q- 1 A N + (-l) N q- 3 
q + q-i{ qA N + (-1) V A N - (-1) N 



a/3 



(5.3.21) 



SajSps (5.3.22) 



-(-1) 



N 



1 

1 



01 



a/3 



where 



A = q 2 + 1 + q- 2 . 

Again the q = 1 limit of this formula recovers the result of H . 

With this expression for the norm, the spin-spin correlation functions defined as 



in the VBS states can be computed by breaking down the numerator into 

(Vs «j . ^ajft )%-7i+i e ft« 3 +i ( n 7^T* ' ) 
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and applying ( 5.3.22| ). We shall display only the result for the special case a = 7 and (3 = 8, where the 
nonvanishing correlation functions are 

< S+SJ >n = (-ly^iq + q- 1 )^-^- 1 - (-iy q - 1 bA N -^ - (-l^-ig&A* -1 

+ (-l) N -i+ l c]/(A N -(-1) N ) 

< S+SJ- >i 2 = (-iy- 1 (? + g - 1 )[oA JV -'- H - 1 - (-1)*?-^*-* + (-l) JV - J ?" 1 6A i - 1 

_ (-i)^+v 2 c]/(A Ar + g - 2 (-l) Ar ) , 

< S+SJ >2i - (-iy- 1 {q + q- 1 )[aA N -i+ i - 1 + {-ifqbA^ - (-l^-iqbA*- 1 

_ {-l)»-i-"q*c]/(A N + q 2 (-l) N ) , 

<S+SJ >22 = {-iy- 1 {q + q- 1 )[aA N -i+ i - x + {-l) i qbA N - i + {-l) N -"q- l bA i - 1 (5.3.23) 

+ (_l)JV-J+i c ]/( A iV_(_ 1 )JVj j 

<S'fS/> aa = (-l)^(g + g- 1 ) 2 (A Ar -^- 1 -(-l) Ar A^)/(A Ar -(-l) JV ), 

< > 12 = (-1) J - 4 (<Z + T 1 ) 2 {q- 2 A N ^+ 1 - 1 + (-1) N A3-*) /{A N + q- 2 {-l) N ) , 

< S*SJ >2i = (-1) J - 4 (<Z + T 1 ) 2 (^A^" 1 + (-1)^'"*) /(A N + q 2 {-l) N ) 

where 

a = g 3 + 2 + <T 3 

6 = q 2 — q + q^ 1 — q~ 2 

c = q + q- l -2. 

Before interpreting these formulae, we first examine the role of the VBS in the spectrum of the hamiltonian. 
For q real, one can extend the proof of the q = 1 case and show that the eigenvalues are always nonnegative, 
and VBS are the only ground states. In the infinite N limit, the ground state is unique with massive 
excitation. Hence, the model is noncritical with spin-spin correlation functions in the VBS given by 

<s+sj>=<sjs+> N =°° (-A)-^v + 2 + 9 -W + i + <r 2 ), (5 3 24) 

<S*SJ> N =°° (-A)-3+*( q 2 + 2 + q - 2 )/( q 2 + l + q - 2 ) . 

The correlation length is therefore l/ln(q 2 + 1 + q ~ 2 ) and notice that the nonisotropy of S+ and Sf in the 
hamiltonian ( 5.3.2) due to the quantum group symmetry is manifested in the above. Only when q — > 1 
where su(2) symmetry is present will isotropy in the spin components be restored. 

Recall that the q — 1 model belongs to the more general antiferromagnetic fluid phase or disorder flat 
phase (DOF) 1 33 . It can likewise be demonstrated that for real q the VBS ground states have the type of 
long range order and disorder associated with DOF phase. The various correlation functions introduced in 
fj3f that distinguish the DOF phase can be calculated. We list, in particular, the density-density correlation 
function 

< {St) 2 {Sjf > VBS N =°°< {S;Y >vbs< {Sjf >vB S Ar = ° (g2 + 1 4 +g _ 2)2 (5.3.25) 
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which confirms that spin positions are completely uncorrelated, and the correlation function which exhibits 
antifcrromagnetic ordering, 



G.U - i) =< Ste^-^) S , >VBS ^oo {q2 + i 4 +q _ 2)2 . (5.3.26) 

The lack of distance dependence shows that AF spin order is perfect. 

For q = e 11 , 7 G R, the configuration space belongs to (C 3 )^ and the reasoning that led to the proof of 
massive excitations for real q no longer holds. The eigenvalues can in fact be negative. 

Numerical check reveals that for finite N the VBS continue to be the only ground state for 7 e [0, 5) and 
7 G (t^ > §]• I n ^ ne fi rs t domain, q 2 + l + q~ 2 > 1, so we expect the properties for real q to be still qualitatively 
valid, with massive excitations and a kind of DOF phase. The second domain has \q 2 + 1 + q~ 2 \ < 1 so the 
behavior is now possibly different from the real q model, in particular, it is not sure whether excitations are 
still massive. For 7 6 (^, there are negative eigenenergies so we certainly expect different properties. 

That another eigenenergy crosses the value zero at 7 = ^ can be shown using U,su(2) symmetry. Indeed 
the projector (Q— 1)(Q — 2)P2 when restricted to type II representations satisfies the Temperley Lieb algebra 
with0 

2 ' 2 — €■% (5 3 27) 

and ef = 2e» . 

The same algebra in the spin-i representation given by 

e, = 2P (5.3.28) 



with 2Po acts on C 2 ® C 2 has "q" parameter of the quantum group given by "q" — 1 or = 2). The 

type II spectrum of 2P 2 at 7 = -| and the entire spectrum of 2P at "q n = 1 share the same set of eigenvalues. 
Moreover the q-dimensions (defined as (2j + l) q ) of the spin sectors from the two representations which share 
the same eigenvalues are equal. In particular, the zero eigenvalue, which occurs in the highest spin sector 
(j = y) of the spin-i representation, has g-dimension 

(2j + l)i - (JV + l)i =N + 1 , 



while in the spin-1 representation, the contribution to the g-dimension of the VBS states, which have j = 0, 1, 
amounts to 

(2-0 + l), + (2-l + l) g = l + g 2 + l + g - 2 =3<7V + l ; for TV > 2 

implying that new zero eigenvalues must emerge for spin chains with N > 2. As an example, at N = 3, we 
find a new zero eigenstate with j ' = 2, the g-dimension of which 

(2-2 + 1), = 1 

adds to the above giving the total contribution 

4(= N+l) . 

The crossing of eigenvalues (and also at j) can also be explained. The j ' = 2 spin representations are 
then type I representations and (Q — 1)(Q — 2)P 2 vanishes when restricted to type II representations. Thus 
all type II eigenvalues vanish. 
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Despite the fact that quantum group symmetry implies additional zero eigenstates have to emerge at 5 
and we only have numerical support that this does not happen outside the domain [5, 

We did not get definite numerical evidence for possible critical properties in the domain Notice 



however the special value 7 = § where the P2 projector line ( 5.3.1 ) meets the TL line, so there we have 
criticality with c = — 2. 

5.4 The T 2 Potts model 

We have discussed the phase diagram of the quantum spin chain because it is the simplest and has most 
immediate applications. It is not always easy to discuss the relation of this study with the two dimensional 
Potts model. Clearly the above hamiltonians, although considered so far as acting on spins, can be rewritten 
as Potts hamiltonians using the appropriate representation of the projectors discussed earlier. It is reasonable 
to hope that the physics of these hamiltonians is the same as the one of a two dimennsional strongly 
anisotropic Potts model whose elementary transfer matrix reads 1 + eH. This correspondence is enough to 
apply to the quantum spin chain duality arguments deduced for a two dimensional (not necessarily isotropic) 
Potts model. However when couplings in two directions take comparable values, it is not clear whether the 
physics will or not be qualitatively different. This is especially true in our case where there are both 
ferromagnetic and antiferromagnetic interactions. In the integrable cases however, one can usually connect 
the physics for different isotropics by changing the spectral parameter, and exact solutions usually show that 
properties are the same provided this parameter runs in a certain range. Let us write the isotropic interactions 
associated with the three integrable lines discussed earlier ( they will be recovered in the hamiltonian limit 
u = 0, e = — 1. The case for e = 1 can similarly be studied. ) 

The Boltzmann weight associated to the fundamental block Qi has physical expression given by 

Wabcd = exp£ afccd (5.4.1) 

where the interaction energy 

£ = K {Q-l)P + lC l {Q-l)P l . 
Written in terms of the four sites a, • • • , d, the energy reads 

Sabcd = (/Co - 4/C0Q- 2 -{lC Q -21C 1 )Q- 1 (5 ab + S cd )+>C 1 Q- 1 (5 ac + 6 bd + 5ad + 5 bc ) (5.4.2) 

+ (/Co — 1Ci)5 a b5cd — fcl(Sabc + 5 a bd + $bcd + ^acd) + K-lQ^abcd , 

which shows that the various interactions; nearest neighbors, next to nearest neighbors etc. can either be 
ferromagnetic or antiferromagnetic depending on the values of the coupling constants /Co and K\. The dual 
model involves similar expressions. 
The TL integrable line is given by 

Q- 1/2 h = o 

which translates into 

, £° = ^ 2 /(Q-D, (5A3) 

and iCx = . 



using ( p.2.8] ) and /q = 1. In this case the energy expression becomes a product of Q —8 a b and Q ~~S cc i with 



/Cq being the overall coefficient. The ferro- and antifcrro-magnetic nature of the interactions therefore depend 
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only on the sign of /Co, which flips at Q = 1. For Q > 1 the model is characterized by antifcrromagnetic 
nearest neighbors interaction and ferromagnetic four sites interaction S a bS c d, and the converse for Q < 1. It 
should be noted that the point Q = 1 corresponds to Q' = of the Ti ( or standard ) Potts model via the 



Temperley Lieb algebra. For the Ti Potts model on the TL integrable line ( 3.6.7), Q' = is precisely the 
point that divides the ferro- and anti- ferromagnetic regimes Q. 
The FZ line is given by 



which is equivalent to 



/C = 
£i = 



Q- 1/2 / 1 = ^ TT , (5.4.4) 
ln(Q + Vg-l)-ln(Vff+l) (5A5) 



where the coupling constants at Q = 1 and Q = 2 are defined by continuity. It is easy to see that both /Cj's 
are nonnegative functions of Q. However /Co becomes complex for Q < (3 — y5)/2 ( or 7 > 27r/5 ) and the 
Potts model beyond that point is not physical. In the domain Q > (3 — v5)/2, the magnetic nature of the 
interaction terms in the energy expression remain unchanged being always ferro- or antiferro-magnetic as 
the coefficients /Ci, /Co — /Ci and /Co — 2/Ci are always positive. 

The same analysis can be performed on the IK integrable line, which is given by 



Q- 1/2 h = -1 + v/4-Q (5.4.6) 

or equivalently 

r _ ln(l + (Q - 1)V4~Q) 

r _ ln(3 - Q + (Q - 2)V4~Q) (5 - 47) 
^ - 7J=2 ■ 

The coupling constants are real for Q G (0.77, 3.80) approximately. As Q increases from 0.77 in this domain, 
the majority of the interaction terms which have coefficient proportional to tC\ changes from ferromagnetic 
to antiferromagnetic or vice versa at Q = 3. In the phase diagram this is the point that divides the various 
regimes of the integrable line. Since /Ci > as in the FZ case, the magnetic natures of the majority of the 
interactions of the IK integrable model for Q < 3 are the same as that of the FZ integrable model. The 
exception being the interactions 5 a bS c d and S a b + S c d, whose respective coefficients /Co — /Ci and /Co — 2/Ci 
change from negative to positive at Q ~ 0.8 and Q ~ 1.9. This similarity in the physical behaviors of 
the interactions supports the conclusion reached in the spin chain study that the two lines are in the same 
universality for Q < 3. 



6 Conclusion 



Tfc Potts models provide a rather different kind of physical models associated with spin-fc/2 representa- 
tions of U g su(2), where the higher symmetry constraints arc encoded in a pattern of complicated interactions 
on a plaqucttc. Besides their "academic" interest we hope they can provide new insight on the physics of 
related solutions of Yang Baxter equation, or universality classes. For instance the standard su(2) symmetric 
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quantum spin chains are related to Q — 4 states Potts models with a mixture of ferromagnetic and antiferro- 
magnetic interactions. The splitting between integer and half integer spin is very naturally observed in this 
picture. A translation invariant quantum spin chain is the anisotropic limit of a four state Potts model based 
on a homogeneous vertex model. For half integer spin (fc odd) this Potts model turns out to be necessarily 
self dual. One therefore expect it, by standard arguments, to be at a critical point. On the other hand for 
integer spin (k even) the Potts model is not self dual, and generically is expected to be in some non critical 
state. This is quite close to the Haldane conjecture. 

The technlogy of quantum groups, Tcmpcrlcy Licb algebras and graphical representations is known under 
other names in the condensed matter literature^]. In particular it was remarked in jlTj that the standard 
Q state Potts model can be related to a quantum spin chain with su(n) symmetry, with the fundamental 
representation on a sub lattice and its conjugate on the other, and n = \fQ. More generally one can speculate 
that systems with quantum group symmetries provide proper analytic continuations of models with ordinary 
symmetries when the rank of the algebra or the size of the representation assume "intermediate" values. For 
instance the spin 1 U g su(2) model, or equivalently the T2 Potts model, can be related to a quantum spin 
chain with su(n) symmetry, once again n = \[Q, but with the adjoint representation on every site. This is 
because (3) g = n 2 — 1. As shown in section 5 of this paper, the phase diagram is rich. In particular several 
critical lines and massless phases are met in the continuation from n = 2 (q = 1) to n = (q = i), which 
is of interest for the quantum Hall effect ][l7| . As the representation gets more complex we expect this 
continuation to "go through" a more and more complicated phase diagram. 

Finally we remark that the quantum group symmetric models are also a particular example of anisotropic 
quantum spin chains. From that point of view, the last paragraph in section 5 represents an extension of 
the valence bond method to a particular anisotropic situation. 

Acknowledgments: I. Affleck, B.Nienhuis and N.Read are thanked for useful discussions. 
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A Boltzmann Factor of Potts model 

The Boltzmann weight of the T 3 Potts model contains interactions between all the six sites a, . . . , /, the 
explicit form is given by 

W(u)aebcfd = (0 — iy 2 {S e f - (Q - l)5 e f(Sbd + Sac) - (8 e fa + S e fb + S e fc + S e fd) + Sabef + S a def 
+ Sbcef + Scdef + Q(S bd ef + S a cef) + (0 — l)(SbdS a ef + S bd S ce f + S ac Sbef + S ac 5 de f) 

- Q(S a bdef + Sacdef + Sabcef + Sbcdef) ~ Q(Q ~ ^){S aC Sbdef + hd&acef) 
+ (Q~ l) 2 Sac5bd5ef + Q 2 Sabcdef} 

+ Q- 1/2 h(Q - 1)- 2 {1 - (S cf + 5 fd + S he + S ae ) - (Q - l)(5 ab + S cd ) + Q(5 abe + S cdf ) 

+ S ae S c f + S ae Sdf + S be 5 cf + SbeSdf + (0 - l)(S be S c d + S ae 5 cd + 6 ab S cf + S ab S d f) 

+ (Q- l) 2 5 ab S cd - Q(S ae S cfd + 5 be S cfd + S cf S abe + S df S abe ) - 0(0 - l)(S ab S cdf + S cd S abe ) 

+ Q 2 SabeS c df} 

+ f(Q - 1)~ 4 {2 - 3Q + (3Q - 2){6 be + 5 cJ + 6 df + S ae ) + (Q - l){S bf + S de + S ce + S af ) 
+ S e f + (Q - l) 2 (S bd + S cd + S bc + 5 ad + S ab + S ac ) - Q(S bef + 5def + S ae f + S cef ) 

- (Q — l) 2 (S c de + Sabf) - 0(0 - l){$bde + S b df + S bc f + S a de + 5bce + $adf + S a cf + 5 ace ) 

+ (1 - Q - Q 2 )(S cfd + S abe ) - (Q - l) 3 (S bcd + S abd + S abc + 5 acd ) + (2 - 3Q)(S be S cf 

+ SbeSdf + S ae Sdf + S ae S c f) - (Q - l)(S de S c f + 5 ae S b f + 5 d fS ec + S be S a f) 

- (Q - l) 2 (S b dS c f + S be S cd + 5 bc 5df + S ae Scd + S bc S ae + S a dS be + S a bSdf + S ad S c f + SacSbe 
+ S ac S d f + SabSfc + S bd S ae ) + Q 2 (S bde f + 5 ace f + S bce f + Scdef + S a bef + S a def) 

+ Q(Q— l) 2 (S bcde + Sacde + Sabdf + S a bcf) + Q 2 (Q - l)(S a bce + S a bde + S bcd f + 5 ac df) 
+ (Q 2 + Q - l)(S c dfSbe + ScdfSae + S a b e Sdf + S a beS c f) + Q(Q - l)(Sbd e S c f + S bd f5 ae 
+ SacfSbe + S ace Sdf + S a d e S c f + S bc fS a e + S bce S d f + S ad fS be ) + (Q — l) 2 {S abe S c d + S cd fS ab ) 
+ S a bcd + {Q ~ 1] ) {Sabc S d f + SacdSbe + SdcbSae + S a bdSfc) 

- Q 3 (Sbcdef + Sabdef + S a bcef + S a cdef) — Q 2 (Q — ^) 2 (Sabcde + S a bcdf) - Q{^Q — l)S a beS c df 
+ (Q - l) 2 (SbdSf c Sae + S bc SfdS ae + S ad S be 8f c + 5 ac 5 be 8f d ) + Q 4 5 ab cdef 

- Q 2 (Q - l){5 a bdeScf + SbcdfSae + SabceSdf + SacdfSbe)} 

+ Q~ 1/2 9(Q - iy 4 {Q 2 + (1 - 2Q)S ef - Q(Q - l) 2 (<5 ac + 5 bd ) - Q 2 (S ae + S be + S cf + Sdf) 

- 0(0 - l)(*ed + Sec + S b f + S af ) + (Q 2 + Q - l)(S efd + S e fb + S e fc + S e fa) + Q 2 (S c fd + S a be) 
+ 0(0 - l)(^6/c + Sbec + Sadf + S a de) + Q 2 (0 - l)(^6ed + S b f d + S ac f + 8 aC e) 

+ (0 - l) 2 (<5 a /fc + S cde + S ac Sef + S bd 5 ef ) + Q(Q - l) 2 (S bd S fc + S bd S ae + S be 5 ac + 5 ac S d f) 

+ Q 2 {S be Sf c + S ae Sfd + 5 a eSf c + 5 be 5 d f) + (Q - l) 3 (<S Q /<5 6d + 5 ec 5 bd + S ac S b f + 5 ac S d e) 
+ 0(0 - l)(S ec Sf d + S a fS be + S ae S b f + S c fSed) + (0 - l) 2 {S b fSec + S a fS e d) 

+ (0 - l) 4 S aC S bd + Q(l - - Q 2 ){Sacef + Sbdef) + 0(2 - 3Q)(5 cde f + S a bef + Sbcef + S a def) 

- Q 2 (Q — l)(S a cdf + Sbcdf + Sabde + S abc e) - 0(0 — I) 2 (Sbcde + S ac de + S a bcf + S a bdf) 

- Q 2 (SbeS c fd + SaeScfd + S a beS c f + SabeSdf) ~ Q 2 {Q ~ l)(SbdeSfc + SacfSbe + S ace Sfd + S bd fS a e) 
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- Q(Q — 1) {SaceSbf + SbdfSec + 5 e f c 5bd + SbdeSaf + SbefSac + S e fdS a c + 5 a f c 5 e d + SaefSbd) 

- Q(Q - l)(S bec Sfd + SadfSbe + SaedSfc + 5 b cf5ae) + Q S 

abe^cdf 

+ Q (2 — 3Q)5 a bcdef 

- Q(Q — l) 3 {Sbde5ac + SbdfSac + 8 ac f5bd + ^ace^bd) + Q 2 (3Q — 2)(8 a bcef + $abdef + Sbcdef 
+ Sacdef) + Q 2 (Q — 1) (Sabcdf + Sabcde + SbdefSac + SacefSbd + SafcSbde + S a ce$bdf) 

- Q(Q - l) 2 {S ac 5 b eSdf + $bdSae$cf) + Q 2 {Q - l)(<W/$>e + Sbcdf^ae + 5 a bde&cf + Sabce^df)} 

B Loop Model Formulation of the Tk model 

We present in this appendix a brief review of the known loop model formulations of the Tk model. 



B.l 

First we recall the standard graphical representation of the Temperley-Lieb algebra |30| [5lJ with the repre- 
sentation space being a set of strands. The generator e, acts on two neighboring strands and produces the 
following configurations 




while the identity leaves the strands unaltered 



1 




With this definition, the algebraic relations ( 2.6 ) are represented as 



and 



which are easily seen to be satisfied in this representation provided every loop is given a weight y/Q. Such 
a reformulation was rediscovered many times, in particular in pT^ ] using valence bond language. 

This gives a geometrical reformulation of the six-vertex model (known also as loop model formulation) 
where the vertex x\l + &2i-i or 1 + are replaced by the graphical combinations 
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+ x 2 




and the lattice is accordingly covered by a collections of closed loops. From the Poots model point of view, 
these are the surrounding polygons of clusters high temperature expansion. 

Recall that there are some subtleties about the models correspondence due to boundary conditions. 

Using the fusion procedure, a loop model formulation can be given to the T k vertex model O, fL5|. 



Graphically, this is done by replacing each of the six-vertex in ( 3.1.4 ) by one of the above configurations. 
Eac 

by 



Each spin- 1 vertex therefore acts on 2k strands, the symmetrizer Sk that acts on k strands is represented 



Sk 



Figure(BO) The composite fc-strand that denotes the symmetrizer Sk 



which is a composite object given by ( 3.1.3 ). As an example for k = 2 



S2=1 -(k e 



has graphical representation 



-1/(2) 




The internal vertices rk(ui), ■ • • , r^u^) that are inserted between the four copies of Sk represented above 

,2 

produce in general 2 K configurations. However, configurations that have any two strands that originate from 
the same symmetrizer joined together have vanishing weight, since this implies the presence of the factor 

PiPv = 

where P\ comes from the symmetrizer and Pq cx e from the fusing of the two strands. Hence for k = 2, the 
nonvanishing configurations are those shown in fig.(Bl). 




Figure(Bl) The nonvanishing strand configurations of the 
loop model formulation of the r 2 model 
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In general, the Tk vertex is replaced by k + 1 nonvanishing strands configurations. (For k — 3, see fig.(B2)) 
These strands are again the surrounding polygons of the fused Potts models introduced earlier. In this fused 
loop model, a closed composite /c-strand obtained by fusing the individual strands in fig. (BO) into loops 
carries a weight (k + l) q . 

In such a formulation the numbers of degrees of freedom at each vertex is vastly reduced, but one has 
instead to deal with nonlocal quantities. 




Figure(B2) Nonvanishing strand configurations for the T3 model 



B.2 



For k = 2, there exist other loop model formulations. The first is due to the observation that] 



(q 2 + q- 2 )Pi(l, i + 1) + q(q 3 - q- 3 )P (i, 1+1), 



(B.2.1) 



e t = (q 2 + l + q- 2 )P (i,i+l) 

satisfy the Birman Wenzl Murakami (BWM) algebra|lj|. The latter contains the Temperley-Lieb algebra 
generated by e$ with = q 2 + 1 + q~ 2 as a subalgebra, and 6,; satisfies the braid group relation 

bib i± ibi = b i± ibibi±i . (B.2. 2) 

They have graphical representations defined by the action on two neighboring strands; 



1 : ZZ 



h: = 



d : 



x 
>< 



Besides the Temperley Lieb and braid group relations, these generators satisfy some other algebraic relations 
which can be represented graphically. Most of these relations are then straightforwardly expressed by regular 
isotopy of the diagrams. The others are : 
1.) The first Rcidcmestcr move 



X>-' > 
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produces a factor g 4 . 

2.) The relation 6, — b^ 1 = (q~ 2 — q 2 ){l — ei) holds, which can be represented graphically as 




-2 




X 



3.) A loop carries a weight (3) 9 . 




In terms of these generators, the vertex given in 



( |3.2.1 ) is written as 



1 + <T 2 /i + (/o + q 2 fi)ei - fibi 



(B.2.3) 



where 



h = h/VQ ■ 



B.3 



The two above mappings have the drawback that they involve dense loop coverings of the lattice. An elegant 
way of mapping the T2 model to a "dilute" loop model is given in M. In a first step one uses the edges of 
the vertices that carry the states |± > to form oriented loops, whose direction is given by the spin arrows, 
while edges with the state |0> are regarded as unoccupied. This gives an oriented dilute loop reformulation. 
The problem then is to find under what circumstances one can get rid of the orientations. The simplest 
way to find a correspondence between an oriented and an unoriented loop model is to suppose that in the 
unoriented model loops have a fugacity. This fugacity can be obtained by a sum of local contributions if one 
gives arbitrary orientations to the loops and sums over all possible orientations, provided a phase factor e ±1 
or e ±1 has been assigned to every turn as follows 



For a lattice which has the geometry of a plane, a closed loop has weight e 2 e 2 (e~ 2 e~ 2 ) if the orientation is 
anticlockwise( clockwise )and so gets fugacity (weight) n equal to e 2 e 2 + e~ 2 e~ 2 . In addition, to every edge 
of the vertex, one can assign, without altering the partition function, a local phase factor as follows 




e 



e 



-1 




a 



b 



a 



-1 



b- 1 
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where the solid dot denotes the center of the vertex. 

In an unoriented model, one has the following local loop configurations 



XXxxxxxXxx 

U>1 U>2 OJ3 LJ4 U>5 U)Q 

where solid (dotted) strand denotes occupied (unoccupied) edge and Wj's, and for the last configuration the 
way the two strands overlap has no significance. If moreover the loops have fugacity n it is then equivalent 
to an oriented loop model, that is to a 19 vertex model, with weights which are products of u>i and a, b, e, e 





K)0,00 




L04 , 


V+0,+0 


= Vb-,o- 






^b+,o+ 


= ^_o,-0 




u>\e~ a b , 


^+-,00 


= ^00,+- 




aft , 


^-+.00 


= K)o,-+ 




u>2ea~ 1 b , 


V++.++ 


= V 




L0 5 + UJ 7 , 


V+-.-+ 


= V-+.+- 




LJ 6 + UJ 7 , 




V+-.+- 




u> 6 e~ 2 a 2 b~ 2 + u> 5 e 2 a 2 b~ : 




V-+.-+ 




w§e 2 aT 2 b 2 + uj 5 e~ 2 a~ 2 b 


Vo-,-o 






L0 3 



(B.3.I) 



2 

V 1- 1- ^(jC W V ~T~ W^G IX V 

V+o,o+ = ^b+,+o 

where Vijja denotes the vertex weight with in- and out- states being ij and kl respectively. 

We thus see that the natural oriented loop model associated with the 19 vertex model is equivalent to an 
unoriented one provided the weights can be parametrized as above. This gives rise to a necessary condition 

y+0,+0 1 V+-,oo (y y \ ^+o.+o 1 y 

VoT^T + v~+w = y y +-<-+ _ v ++:++> v^~o+ + v +-+- (B 3 2) 

Vb+,o+ , V- + ,oo (y y \ Vb+,o+ it/ \ • • I 

( v +--+ - v ++,++) v^; + 

This holds in particular for the Y2 vertex model (ie when the 19 vertex model has U q su(2), for which the 
correspondence between the parameters Wj, a, 6, e, e and /o, /1 is given by 

Wl = (l + g-W^a + g 2 /!) 1 / 2 , 
^2 = (/o_+9- 2 /i) 1/2 (/o+<Z 2 /i) 1/2 , 

^3 = -/l , 

^4 = l + /o + (Q-3)/i, 

w = (l + g 2 /i)(l + g- 2 /i)(/o + /i) 
~ /o + /1 + 7oA + W- 3)/ 2 ' 

(/o + g 2 /i)(/o + g- 2 /i)(i + /i) 



(B.3.3) 



W6 
UJ-J 



fo + A + /0/1 + (Q- 3)/ 2 
A 2 + (Q - 3)/o/i + Ml + g 
h + A + AA + (Q - 3)A 2 



and the loop fugacity reads 



rj _ fj2 (l + g- 2 /l)(/0 + g- 2 /l) ,„-2 (l + g 2 A)(/0 + g 2 /l) , R „,, 

9 (i + g 2 A)(/o + 9 2 A) 9 (i + g- 2 A)(/o + <z- 2 A) ' 
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In this loop reformulation there are more degrees of freedom at each vertex than in the first one we discussed. 
Not all edges are occupied. The fugacity depends on fo and f\. 
Consider now the case 

/i = -l, / = 0. (B.3.5) 
The nonvanishing weights after rescaling become 

^2 = W 3 = U> 5 — 1 , 

and uj\ = L04 = —{q — q 1 ) 

the vanishing of oj§ and and LU7 implies that, if vertices are "expanded" as 



(B.3.6) 





such that the entire lattice becomes honeycomb , every edge can at most be occupied by one strand. Moreover, 
since the nonvanishing weights satisfy the relations given above, the model belongs to a subset of the class 
of loop model where loops do not intersect and the only parameters are loop and vacant site fugacities]55| 
{■n,u>i}. In this case, the two parameters are related as 

n = 2- (2-uj) 2 . (B.3.7) 

The integrable lines IK, TL and FZ play distinctive roles here too, they arise as a result of the restriction 

IK, TL : uj-j = 

FZ : vertex weights being invariant under reversal of all arrows and n=2. 

For general fo, fx, the loops can be interpreted as high temperature expansion of an O(n) model which 
has Boltzman weight 

LUi + Ul(Si ■ S k + Sj ■ Si) + LU 2 (Si ■ Sj + S k ■ Si + U 3 (Si ■ Si + Sj ■ s k ) 

w 5 (s*i • Sj)(s k ■ si) + uj 6 (si ■ s k )(sj ■ si) + u> 7 (si ■ si)(sj ■ s k ) 
where Si 's are n-component vectors situated on the edges of each vertex and are normalized 
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